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AN APPROXIMATE SOLUTION OF THE EQUATIONS OF MOTION 
FOR ARBITRARY ROTATING SPACECRAFT 

By Peter Ralph Kurzhals 

ABSTRACT 

The determination of the motion of rotating spacecraft, such as 
manned space stations and spinning satellites, requires the solution 
of the spacecraft’s equations of motion with varying disturbance torques 
and mass distributions. The numerical integration of these equations 
on high-speed computing equipment can give only limited information on 
the effects of disturbance and spacecraft characteristics, and cannot 
provide the physical insight needed for an analysis of the spacecraft 
motion. An approximate solution of the governing equations which would 
yield a direct assessment of the effects of applied disturbances and 
would lead to a clear understanding of the motion mechanics could thus 
be particularly useful. 

This dissertation comprises the development and application of an 
approximate analytical solution for the motion of arbitrary rotating 
spacecraft with variable disturbance functions. The solution is based 
on the assumptions of small changes in the spacecraft inertia character- 
istics, body rates, and Euler angles. The rate and attitude errors, 
resulting from the application of disturbance torques, are described by 
complex pseudovectors and the governing spacecraft equations are reduced 
to linear differential equations in terms of these error vectors. Solu- 
tions are obtained for the steady spinning mode and for a spinup and 
despin mode. 
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on high-speed computing equipment can give only limited information on 
the effects of disturbance and spacecraft characteristics, and cannot 
provide the physical insight needed for an analysis of the spacecraft 
motion. An approximate solution of the governing equations which would 
yield a direct assessment of the effects of applied disturbances and 
would lead to a clear understanding of the motion mechanics could thus 
be particularly useful. 

This dissertation comprises the development and application of an 
approximate analytical solution for the motion of arbitrary rotating 
spacecraft with variable disturbance functions. The solution is based 
on the assumptions of small changes in the spacecraft inertia character- 
istics, body rates, and Euler angles. The rate and attitude errors, 
resulting from the application of disturbance torques, are described by 
complex pseudovectors and the governing spacecraft equations are reduced 
to linear differential equations in terms of these error vectors. Solu- 
tions are obtained for the steady spinning mode and for a spinup and 
despin mode. 
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The solutions for the spinning mode consider the effects of initial 
errors, external torques, and instantaneous and periodic mass motions 
within the spacecraft. The resultant errors are presented as error 
component time histories and as traces of the complex error vectors. 

Upper bounds of the error magnitudes are deduced from the error vectors. 
Both the general case of nonsymmetric spacecraft and the special case 
of spacecraft rotating about an axis of symmetry are examined. 

Periodic mass motions within the spacecraft are shown to have 
significant effects on the spacecraft motions and can produce errors 
several times greater than the errors predicted for "worst-case" 
instantaneous mass motions. Instability trends of the errors are also 
found when the spin axis becomes an intermediate axis of inertia during 
a mass motion and when the motions occur at the precession frequency A. 

The effectiveness of several control techniques is investigated for 
the approximate governing equations. Pure rate control and rate plus 
rate integral control are found to be acceptable for damping of the 
rate and attitude errors produced by mass motions and other internal 
disturbances. Rate plus attitude control is, however, needed for the 
eli m ination of possible residual attitude errors due to external 
disturbances and for the reorientation of the spacecraft. The implemen- 
tation of the control techniques is discussed for reaction wheel, control 
moment gyro, and reaction jet systems. Actuator commands and the required 
control system weights are developed. 

A comparison of the analytical solution and the exact solution 
obtained from numerical integration of the complete equations of motion 
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is used to establish the adequacy of the approximate solution. The 
applications of the analytical solution are illustrated for a manned 
orbital research laboratory and a large spinning space station. 

The solutions for the spinup and despin mode are employed in the 
optimization of spinup and extension techniques for cable- or strut- 
connected spacecraft modules. Fuel savings of about 22 lb per spinup 
and despin cycle of the manned orbital research laboratory can be 
obtained by a continuous-thrust extension. 

The analytical solution shows that a simple and valid interpreta- 
tion of the spacecraft motions is possible for a large number of applied 


disturbances. 
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IV. SIMBOIS 


Aj complex coefficients of forcing function, (B-l) and (30) 

a power conversion factor, (254) 

a,b,c,d,e,f characteristic coordinate for the error traces 

complex coefficients of total rate error relation, (43) 
Cj complex coefficients of total angular error 

relation, (44) 

B2>D 2 ,Dj,D^,I^ spinup parameters, (295) and (297) 

E 1 ,E 2 ,E 3 ,E4 forcing function coefficients for variable inertia 

products, table 4 and (112) 

Eg,Ey solution function coefficients for variable inertia 

products, (ll3) and (115) 

F complex forcing function, (30) 

F complex solution function for the body rate errors 

F complex solution function for the Euler angle errors 

fx>fy forcing functions for the controlled spacecraft, 

(213) and (214) 


G complex actuator torque 

G s actuator stall torque 

g complex control torque, (209) 

H rigid-body angular momentum 

I moment or product of inertia 

I 0 spacecraft moment of inertia without moving particles 

I rz radial inertia product, I X2 + iIy Z 
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K jx' K jy 

k jx> k jy 


specific impulse 
total impulse, (28U) 

control gains for the X and Y axes, (209) 
nondimensional control gains for the X and Y axes, (210) 
nondimensional external moment, (A-30) 
distance between mass centers of manned and counter- 
weight modules 

effective spinup moment arm, ( 29*4-) 
external moment, + iMy 
spinup moment 
mass of moving particle 



frequency for moving particle 

constant defining center of mass change, (A-25) 
position vector of moving particle measured with respect 
to the origin of the fixed coordinate system 
position vector of origin of spacecraft coordinate 
system measured relative to the origin of the fixed 
coordinate system 
damping ratio 

position vector of moving particle measured with respect 
to the spacecraft coordinate system 
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s 

sgn 

T 

T z 

t 

t D 

U(t) 

V 

W 

X,Y,Z 

x,y,z 

u,v,w 

a 

a g 

a I 

5(t) 


position vector of coinposite mass center measured with 
respect to the spacecraft coordinate system 
Laplace transform variable 

signum function, denoting the sign of the characteristic 
coordinate or variable 
constant step or impulse torque, T x + iT y 
cross coupling torque applied to spacecraft by control 
system 
time 

time constant 

unit step function, (B-4) 

effective forcing function 

weight 

reference axes 

scalar components of "r along the X, Y, Z axes 
nondimensional scalar components of *r along the 
X, Y, Z axes (a-30) 

complex angular position error, cp + iG, figure 4 

limiting gyro gimbal angle, (272) 

complex inertial position error, (l4) 

argument of Bj 

argument of Cj 

impulse function, (B-5) 
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e,cp,t 

A,A k 


M. 

P 

a 

T 

$ 

n 


^xy 

con 

“fc 

COx,^,^ 


Subscripts: 

c^ d.^ G 

c 

CM 


quantity denoting the magnitude order of co^, d>p, 0 ^, 

2 

0 , Lp, pp, and the nondimensional inertia terms 
for the moving particles 

nondimensional moment or product of inertia, (A-29) 
modified Euler angles, figure 2 

precession rate parameters, (24), positive when I z 
is maximum inertia and negative when I z is minimum 
inertia 

nondimensional mass, (A-29) 
nondimensional force* (A-30) 
constant positive spin rate, ( 15 ) 
nondimensional time, (A-29) 

angular coordinate used in total error traces 

total angular rate vector of spacecraft axis system 

— > 

scalar components of ft along the X, Y, Z axes 
complex rate error, % + ifty, figure 3 
damped natural frequency 
synchronous wheel speed 

— > 

nondimensional scalar components of ft along the 
X, Y, Z axes 

value for corresponding characteristic coordinate 
counterweight module 
extension at constant momentum 


CR 


extension with constant spin rate 


CT 


d 

F 
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G 
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i 

j 

k ,1 
lim 
M 
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max 
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o 

P 

P,<1 

P<1 

SU 

s 

R 

r 
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extension, with continuous thrust 

disturbance 

fixed coordinates 

final value after spinup 

gyro 

gyro gimbal 

intermediate coordinates 

intermediate value before extension 

value for jth term or mass where j = 1, 2, 3> • 

component for X or Y axis with k / 2 

upper bound 

momentum 

manned module 

maximum value 

summed value for moving particles 

value referred to origin of spacecraft body axes 

initial value 

power 

component for X, Y, or Z axis 
component for XY, XZ, or YZ plane 
spinup fuel 

value for spacecraft mass center 
reaction control 
residual value 



- 15 - 


T total 

x,y,z component for X, Y, or Z axis 

xy,xz,yz component for XY, XZ, or YZ plane 

W reaction wheel 

A dot over a symbol denotes the derivative with respect to time. 

An arrow (-») over a symbol denotes a vector. 

A tilde W over a symbol denotes the general solution function 
corresponding to initial rate and attitude errors. 

A single (') and double (") apostrophe denotes particular elements 
of a vector component along the spacecraft axis. 

The quadrant for the angles corresponding to the inverse trigono- 
metric functions tan~^( ) is determined by the sign of the numerator 
and denominator of the term in the brackets. When both numerator and 
denominator are positive, the angle is in the first quadrant; when the 
numerator is positive and the denominator is negative, the angle falls 
in the second quadrant; when both numerator and den omina tor are negative, 
the angle falls in the third quadrant; and when the numerator is negative 
and the denominator is positive, the angle falls in the fourth quadrant. 

All square root terms in this analysis are principal, positive 
values. These values may be positive real or positive imaginary numbers. 
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V. SUMMARY 

The assumption of small changes in the inertia parameters has been 
used to derive approximate rotational equations of motion for arbitrary- 
spinning spacecraft in the small angle and rate regime. Complex repre- 
sentations are introduced to define the rate and attitude errors produced 
by applied disturbances, and analytic solutions are obtained for the 
steady spinning mode and for the spinup and despin mode. 

Solutions for the steady spinning mode consider both the uncontrolled 
and the controlled spacecraft motion for characteristic disturbances. 

These disturbances include initial errors, externally applied torques, 
and instantaneous and periodic mass motions within the spacecraft. The 
errors induced by the disturbances are described by the error component 
time histories, and by vector traces of the complex error representations. 
Upper bounds of the errors are developed for the uncontrolled case, and 
the required control techniques and control systems are examined for the 
controlled case . 

Solutions for the spinup and despin mode consider extensible space- 
craft modules connected by struts or cables. Fuel consumption relations 
are derived for several extension techniques, and optimization of the 
extension techniques is shown to yield appreciable fuel savings. 

Comparisons of the analytical solutions and exact solutions obtained 
by numerical integration of the complete equations of motion are found to 
be in excellent agreement, and the applications of the approximate solu- 
tion are illustrated for a manned orbital research laboratory and a large 
spinning space station. 
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VI, INTRODUCTION 

Proposed spacecraft, such as the manned orbital laboratory (ref. l) 
and manned interplanetary vehicles (ref. 2), may use rotation about a 
maximum axis of inertia to provide spin stabilization and to produce 
artificial gravity for the crew. These spacecraft will be subjected to 
variable torques arising from both internal and external sources (ref. 3 ) 
and will undergo wobbling motions as a result of these torques. Since 
the wobbling motions (ref. 4) produce attitude errors (which may affect 
the spacecraft’s power system and experiments) and oscillatory rates 
(which may lead to discomfort and nausea of the crew), an analysis is 
required to determine the magnitude of any such attitude errors and body 
rates for the spacecraft under consideration. 

In order to carry out this analysis, the spacecraft's equations of 
motion with varying inertias and torques must be integrated to define 
the spacecraft response for the anticipated applied disturbances. In 
the past, such a solution has required high-speed computing equipment 
for the numerical integration of the equations of motion and has con- 
sumed a large amount of computer time to assess the effects of a range 
of disturbances for a particular vehicle configuration. 

Because of the rather limited application of these results, an 
approximate analytical solution of the spacecraft’s equations of motion 
would be of considerable value. The closed form solution could be used 
to determine attitude errors and body rates introduced by "worst case" 
type of disturbances and would define instability trends that might 
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result from applied torques. In addition, such a solution would allow 
a direct evaluation of the effects of changes in both the spacecraft's 
configuration and the disturbances on the spacecraft's motion. 

Approximate analytical solutions of the equations of motion for an 
arbitrary rotating spacecraft may be obtained for linearized governing 
equations and Euler angle transformations. A number of such solutions 
have been obtained for the simplified equations of motion corresponding 
to symmetric or near-symmetric spacecraft. Leon (ref. 5) and Thomson 
(ref. 6) have developed attitude and rate relations for spinning near- 
symmetrical bodies by considering a vectorial representation of the total 
errors. Thomson and Fung (ref. 7) have also investigated the stability 
of near- symmetric spinning space stations and have defined regions of 
instability for an example vehicle. In addition, Hackler (ref. 8), 

Buglia (ref. 9)> and Loebel (ref. 10) have derived expressions for the 
attitude and rate histories of symmetrical spacecraft by linearizing the 
equations of motion. 

Several analytical solutions for a nonsymmetric spinning body with 
constant inertias have also been obtained. Exact solutions for a torque- 
free body were developed by Routh (ref. ll) and MacMillan (ref. 12) in 
terms of Poinsot's construction and elliptic functions and by Whitbeck 
(ref. 13) in terms of a phase plane approach. An approximate method 
which shows good agreement between the nonlinear and linearized results 
for a vehicle under applied torques was presented by Suddath (ref. 1^). 

Various other analytical approximations are discussed in the liter- 
ature (refs. 15-20). Existing solutions, however, have considered either 
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very special cases of nonsymmetric spacecraft or have been restricted, to 
particular symmetric or near- symmetric spacecraft with specified, distur- 
bances, Such results cannot be applied to the general case of a nonsym- 
metric spacecraft with varying products of inertia and applied torques, 
and offer little information on the properties of the motion of such 
spacecraft. Furthermore, the form of these solutions has made the 
determination of upper limits, for the total attitude and rate errors 
difficult since the amount of computational time required to define the 
error boundaries is in general prohibitive. 

The present analysis develops a solution technique for arbitrary 
rotating spacecraft with variable disturbance functions. The complete 
equations of motion for nonsymmetric vehicles are linearized and solved 
with time varying forcing functions and products of inertia. General 
and particular solution functions are determined and are used to generate 
rate and attitude expressions corresponding to the variable forcing 
functions. A complex vector representation is introduced to define both 
error time histories in component form and the total angular and rate 
errors . 

A number of disturbances are considered for both nonsymmetric and 
symmetric spacecraft; and the corresponding solutions are examined for 
the uncontrolled and controlled cases. Upper bounds of the total errors 
are defined and body-fixed and inertial traces of the total errors are 
analyzed. A method of selecting control commands is also presented. 
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VII. PROBLEM FORMULATION 
A, Spacecraft Motion 

The rotating spacecraft will he related to the reference system 
shewn in figure 1. A set of X Y Z axes fixed to the spacecraft is used 
to describe the rotational motion of the spacecraft with respect to a 
set of intermediate Xj Yj Zj coordinates. The intermediate coordinates 
translate without rotation in inertial space, but always remain parallel 
to a set of Xp Yp Zj> axes fixed in inertial space. 



Figure 1.- Reference system for rotating spacecraft. 

The inertial attitude of the spacecraft may be defined by means of 
three modified Euler angles which determine the relative motion between 
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the X Y Z and Xj Yj Zj axis systems. These modified Euler angles, as 
illustrated in figure 2, results frcm three consecutive rotations. The 
first rotation, about the Zj axis, carries the Xj and Yj axes through an 
angle t measured in a horizontal plane. The second rotation, about the 
new Yj axis, then takes the Xj and Zj sixes through an angle 9 measured 
in a vertical plane. Finally, the third rotation, about the new Xj axis, 
carries the Yj and Zj axes through the angle cp measured in an inclined 
plane to give the X, Y, and Z axes. 



Figure 2.- Vector transformation between spacecraft axes 
and intermediate reference system. 


The modified Euler angles can be determined by expressing the 

• • • 

rotations ty, 9, and cp in terms of the angular rates Qy, 
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and H z about the vehicle axes, The vehicle angular rates then can be 
found from a solution of the vehicle moment and force equations. The 
resultant expressions for A*, %> and flz are substituted into the 

Euler angle transformations, which now reduce to differential equations 
in \|r, 9, cp, and t. The solutions of these equations give the attitude 

of the spacecraft relative to the intermediate axes and thus determine 
the angular motion of the spacecraft. 

B. Assumptions 

To make the general nonlinear equations of motion amenable to 
analytical treatment, the assumptions 

°z A 

« % 

sin 9 = tan 9=0 cos 0=1 

sin cp = tan cp = cp cos cp = l J 

were introduced in the moment, force, and Euler angle relations developed 
in appendix A. The further assumption that the nondimensional inertia 
terms associated with any mass particles moving with respect to the 
spacecraft were small was also made to linearize the equations. The 
resultant method of reduction to linear form and the range and validity 
of these assumptions are discussed in appendix A. 

C. Governing Equations 

With the assumptions of the preceding section, the equations of 



motion reduce to 
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The inertia terms are 


ll 

X x = + y m j(yjo 2 + z jo 2 ) - m s(y S o 2 + z so 2 ) 




J-l 

n 


I y = Iyo + y m j( x jo 2 + z jo 2 ) - m s( x so 2 + z so 2 ) 
n 

*z = Izo + y m j( x / + y J 2) “ m s( x s 2 + y s 2 ) 
j-l 


' ( 3 ) 
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I xz = j mjXjZj - m s x s z s 


x yz = ) , m jyj 2 j - "Ws 
J=1 


and, consistent with our assumptions, the moment s of inertia are taken 
as constant in (2) and (3), hut are allowed to vary in ( 4 ), The asso- 
ciated inertia derivatives then become 


n 

x z = 2 7 + yjyj) * m s(x s * B . + y s y 5 ) 


11 

ixz - y + Vj> - m s ( x s z s + z s*s> 


n 

iyz * / + z jyj) - m s(y s z s + z s y E ) 

M J 

where Xj, yj, Zj denote the position coordinates of the mass 
moving with respect to the spacecraft, and 


, y a 

L m s 




„ y a 

U m s 


n 

V 

L % 

J=i 


denote the position coordinates of the composite mass center for the 
spacecraft and the moving masses. 

The spin rate Q z is obtained by integrating (4) as 


n z ~ —-(IzO^zo + J Mz dt + J 7 - xjyj) 

z l M 


- m s (ysXs - x bVb) 


dt 


(8) 


where the first term in the brackets represents the system's initial 
angular momentum and the remaining terms account for changes in the spin 
rate due to the applied torque Mz and to the accelerations of the moving 
masses. 

Solutions to the spacecraft equations of motion may be obtained by 
first determining ftj. from (8) and then integrating (2) and (3) simul- 
taneously to find fljj and fty. 

These solutions can be substituted into the linearized Euler angle 
transformations 


q> = % + GzG 


(9) 


0 « «y - 
♦ = «z 


( 10 ) 

( 11 ) 


from which we have 

♦ 


= J J-<Iz0^o + J Mj dt + J 


~n 

i ’ 


- x jyj) 
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The solution of (2), (3)> (9)> and (10) then defines the motion of the 
spacecraft in terms of the time histories of the body rates and Euler 
angles. 


D. Total Errors 

The body rates ^ and fly are the undesired rate components pro- 
duced by the applied disturbances and will be referred to as the rate 
error components. Similarly, the Euler angles 9 and 0 describe the 
unwanted attitude deviations that result from the application of the 
disturbances. These Euler angles will be referred to as the attitude 
error components. The solutions for both rate and attitude error com- 
ponents follow directly from the preceding section, and are found as 
time dependent components along the body and inertial axes. 

In practice, one is primarily concerned with the total errors. For 
example, the time variation and maximum value of the total angular 
velocity error in body-fixed coordinates must be known to assess possible 
crew discomfort d ue to wobbling motions. The time yariation and maximum 
value of the total angular position error with respect to inertial space 
is needed to determine possible effects on the spacecraft experiments 
and power system. The effects of removal of a disturbance on the 
residual spacecraft motion are also of interest. 

Both the total angular position and the total body rate errors may 
be developed by using a complex vector representation (ref. 5)* The 
total angular rate error fl^y can be obtained by vector addition of 
the body rates fl* and fly, as shown in figure 3. 
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Figure 3*- Vectorial representation of total angular 

rate error. 


Mat hemat i cally , 


ftjjy may be written as 

I g p itan ” 1 

^xy " + ~ \j®x + % 6 


^x 


( 13 ) 


Similarly, the total attitude error a in body-fixed coordinates can 
be considered as the vector sum of the small Euler angles cp and 8, 
as illustrated in figure 4. To transform this pseudovector to the 
intermediate coordinate system, one must rotate the body coordinate 
system through the angle The total inertial error aj is then 


a x . aa 1 * = (<p + lS)e W . 



O'O 
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from (9) and (lO), one may use (ll) to develop the relation 


Oj = (a + if^aje 
= fi* ye 


ifi z t 


i9 z t 


(1Y) 


The magnitude of the rate of change of the inertial attitude error is 
thus equal to the magnitude of the rate error for the small angle 

regime. 

By integrating (lT), one arrives at 


a I 



dT 


(18) 


us the solution for the inertial attitude error vector. The attitude 
error in body coordinates becomes 



and both the attitude errors a,j and a can be directly developed 
from the rate error expression. 

If only total error vectors are desired, the time solution of the 
equations of motion for may be followed by application of (18) 

and (19) to yield a.j and a. If the rate and Euler angle components 
are of interest, the direct solution of the linear differential equa- 
tion (3), (9)> and (lo) is preferable. 

In the present analysis, solutions were first developed in the 
i.' m-iu of time histories for the error components. The transition to the 
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total error form was then made by substitution of the vector expression 
for the resultant rate error in (l 8 ) a nd ( 19 ). 

E. Solution Approach 

The solutions of the equations of motion for arbitrary rotating 
spacecraft can in general be divided into two types, namely those 
associated with the spinup and despin modes and those associated with 
the steady spinning mode. 

The spinup and despin modes may involve the extension and retrac- 
tion of cable-connected counterweight modules and thus could have major 
and rapid changes in the moments of inertia for the spacecraft. During 
these modes other disturbances, such as crew motions and applied torques, 
will necessarily be restricted and only the solution for the spin rate 
and angle, as given by ( 8 ) and (12) need be considered. The efficiency 
of various spinup and despin methods using constant spin rate, constant 
cable tension, or similar schemes can be readily evaluated from these 
equations. 

For the steady spinning mode, the variations of all total moments 
of inertia due to the moving masses associated with a particular crew 
motion are small in comparison with the constant spacecraft inertias 
I X 0 > IyO and IzO The assumption that the total moments of inertia 
I x , I y , and I z retain their initial values throughout the crew 
motion (see appendix A) will be made for this mode. 

In addition, disturbance moments due to crew motions and applied 
torques will now act primarily about the spacecraft X and Y axes. Any 
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torques nbout the spacecraft Z axis can be neglected during a particular 
• I i stu rbance since the resultant change (refs. 8, l 4 , 19) in the spin 
I'oi.e will be small in comparison with the initial spin rate. 

In accordance with these assumptions, one may approximate the spin 
rate by its constant value at the initiation of a particular disturbance 

% - s „ (20) 

for the evaluation of the effects of that disturbance on the spacecraft 
motion in the steady spinning mode. The value a can and will be 
taken as positive without loss of generality. 

Since the spinning mode occurs for the major portion of the space- 
craft lifetime, this mode will first be analyzed in considerable detail 
and several spinup and despin techniques will then be considered in a 
later chapter. 
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VIII. ANALYSIS OF SPINNING MODE 

For the spinning mode, the moments of inertia take on their initial 
values immediately after initiation of the disturbance and remain con- 
stant for the duration of the disturbance. The inertias may thus be 
computed from (5) as 



and the governing equations may now be developed directly from (2) 
and (3). 


Equations of Motion 


33 



and 
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For simplicity of notation, introduce the precession rate parameters 






^ 2 - w 


so that (22) and (23) became 
" 4 c + A % = /^[ x - AyMy + c|ixz " iyz( 


a + Ay) - IxZ^Ay 


n r 

+ £ mj (<r + Ay)(zjXj - x j z j) - a V z jyj - y j z j> 

A 

+ ^jy'j + z j y ,i - yj®j - yj^j) - m sjy ff + V (z s x s - 


” a A z sy s " y E z s) + (z s y s + z sys - ys z s ~ ys z s) 


Qy. + A Cly = + + ^|^-yz + ^-xz^ 11 + ^x) " Xy Z crA x 

v r 

+ ) m s (a + A*)(*jy 3 - yjz'j) + o\(zfa - x^) 

A L 

+ + k A - z j x j - Z A)] " m s[<> + A x)( z s y s 


-V / • • \ / • •• *«• • *• V 

+ a\{z s x G ~ x s z s ) + (x s z B + x s z s *- ZgXg •* z s x s) 


- y s z s) 
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or 


°x + * Ox “ F > 


+ * <V - F y 


(25) 

(26) 


where 


£ " ^y^y + 0 [*xz " iyz( a + \) “ 


n 

+ V m 


+ vtzjxj - x 3 z j) - <r*y( z jyj - yj z j) 

(zjyj + Z.yj - yjZj - yjZj)J - m^a + *y)( 


Z S^£ " X S 2 S ) 


- ^(z s y s - y s z s ) + (* B ys + z s y s - y s z s - y s z s) \ 

J J 

= i~ - W + ^ {yj + < 2cT + - 0 < a + 2 Vyj 

- ^ 2 Vj) + z j (yj + 2o *j - Ai) - *j (yj + Vo) 


- Vj 


- m r 


z s {jfe + ( 2cr + ^y)*s - ^(<7 + 2Ay)y s - 0 2 7yx s ^) 


+ z s{y s + 2ox s - c^yg} - z s {y s + A y x s } - ZgyJ 


(27) 
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and 


F y - — + \Mx + a ji*yz + ixz( 0 + \) m ^yz 0 ^ 

+ V m A k<x + \)(z^ - yjij) + - x^) 

J»1 L 

+ W*3 + V'j - z ^j - z j^)j - m s[( ff + \)( z s y s - y s z s) 

+ a7s ac (z s x s - x s z s ) + (x s z s + x s z s - ZgX s z s x s ) 

= ^:^y + ^ ” Z " (2a + ^)yj - CT (° + 2 \)*j 

+ a*h&3) + z j (xj - 2oyj - cr^xj) - zj (ij - A x yj} - z>jj 


+ m c 


*s{*s - ( 2cr + \)vs ” <*(<* + 2 ^x)x s + o 2 \y s ') 


"*■ z s (*s " 2cry s - c^x s ^ - z s ^x s - A x y s ^ - z s x £ 


Adding ( 25 ) and (26) in quadrature and referring to (13) yields 


(28) 


Ojjy + A fijjy - F 


( 29 ) 


with 


F = F + iF 
x y 


( 30 ) 
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The solution of ( 29 ) is 

&Xy = sin + flxy 0 COS ^t + F (51) 

where F is obtained by replacing functions of t in F by the cor- 
responding particular solution functions given in appendix B. 

The initial conditions at t = 0 are 

^xyo = ^xo + ia yo 

^xy° * + 



The particular contributions of an applied disturbance to the initial 
errors are included in the Laplace formulation of the solution terms. 

Substitution of (32) into (3l) then results in expressions for the 
total error ftjjy and its components and fly. The spin rate fi z 

is found from 

n _ ^zO^zo = 

— ' 11 1 — 0 


and all of the body rates have thus been defined. 

The Euler angle differential equations can be written as 


C|5 + (J 2 *? = fljj + (Xfiy 


0 4* (T 0 ~ Q/y » 


(33) 


and after adding in quadrature 


a + c a = flxy - itf^xy 


(3M 
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Substitution for the rate error yields 

sin At + ftxyo - i(,fl xyo| cos At + F 


a + sin At + [v - 


ioF 

( 35 ) 


which has the solution 



where F and F are obtained by replacing functions of t in F by 
the corresponding particular solution functions of appendix B. 

Initial conditions at t = 0 are 


a D = <p 0 + i© c 


®o “ ^xyo “ iaa c 


and the Euler angle i If can be determined from the relation 


( 37 ) 


t = crt 


( 38 ) 


This completes the development of the Euler angles. 

Since the terms involving the initial errors will have the same 
form for all disturbances, introduce 

'w /v m 4*xVO 

53 sin Xt* + rij^y-Q cos Xt 
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or in component form 



- hi - 

All coefficients in these equations may be evaluated from the initial 
conditions 0xo> Sly 0 , a , cp 0 , and the initial moments of inertia. 

The terms involving the applied disturbances can be similarly put 
into component form, so that the body rates and Euler angles can be 
found by equating real and imaginary parts in 

fijcy = fl X + ifly = flxy + f* 

and (42) 

a = cp + 19 = a + F - iaF 

where F, F, and F are taken from appendix B. 

B. Total Errors 

The total rate and attitude errors may be put into a somewhat 
singpler form by expressing all trigonometric terms occurring in these 
errors in exponential form. Thus, one is able to obtain 

v - X f. B J thelPjt 

0=1 h=0 

where j and h vary over a finite range of integers. The complex 
constants Bj and the real constants (3 j must be evaluated for a 
particular disturbance. 

The total angular error in inertial space, as defined by equa- 
tion (l4) can be similarly expressed as 

aj = (cp + iele 1 * = y J 1 OjtV 7 -^ (44) 

j=l h*) 
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where j and h again remain finite integers with Cj and y ^ 
determined for a specified disturbance. 

In engineering applications, one is also interested in the maximum 
magnitudes of these errors. Since the exact solution for the maximum 
error magnitudes requires an iterative determination of the zeros of the 
magnitude derivative - and this requires considerable computing time - 
an alternate method of defining upper bounds for the errors is preferable 
from the practical standpoint. It is noted from (43) that 


i - 


M ■ 


t g 

[7 

j=l h=0 


y y ^ 


f g 

< V V U 


^ L 

J=1 h=0 


i^t 


* 1 1 NK 

j —1 h=0 


hi = 


- j^xyj 


lim 


and similarly from (44) that 


(45) 


a I - 


= a = 


u v 


y y 


j=l h=0 
u v 


s I I |c/| 

,3=1 h=0 


u v 


<- 1 I hi 


t h | = la 


lim 


3=1 h=0 


(46) 
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These upper bounds provide limiting values of the error magnitudes which 
are adequate for assessing the effects of particular disturbances. More 
accurate estimates of the absolute maximum errors and their directions 
can be obtained from polar plots of the complex errors if desired. 

C. Characteristic Disturbances for 
Nonsymmetric Spacecraft 

Most disturbances acting on rotating spacecraft may be approximated 
by impulsive torques, step torques, step products of inertia, or variable 
products of inertia. For example, docking impacts and attitude control 
moments can be represented by impulsive or step torques, while crew or 
cargo motions would result in either step or variable products of 
inertia. Other externally applied torques (such as the sinusoidal 
gravity gradient moments) are dependent on the particular spacecraft 
and orbital characteristics and cannot be defined without selecting a 
specific vehicle and orientation. 

The effects of characteristic disturbing functions on the spacecraft 
motion are presented in this section. Time solutions for the Euler angle 
and body rate errors are developed for arbitrary constant moments of 
inertia, and upper bounds for these variables are given. 

1. Impulsive Torques 
a. Time Histories 

Docking impulses caused by resupply and rendezvous vehicles or 
micrameteorite hits may result in impulsive torques acting on the space- 
craft. These torques can be written as 


186 
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M a M x + iMy = (T x + iTy)&(t) 


and the corresponding forcing function is 


: |i x 5(t) - \TyB(t) + ijT y 5(t) + \T x 6(t)J 


For arbitrary initial conditions, the total rate error may be found 
from (27), (42), and (48) by using the solution functions given in 


table I. The results are 


^xy = °xy + 


i[ T: 


T A 

' x cos At =■—=- sin 

* A 


ln xt ] + if 


VSc 

Ty cos At + — ^ — sin At 


where ft^y is given by (39) and (4o) . 


The attitude error can be similarly determined as 


a = a + 




.(cos At - cos crt) + T_(~ sin At + sin ort^ 


- i T 

L. 


T x (cos At - cos crt) - sin At + sin crt 


with a determined from (39) and (4l), 


b . Total Errors 


Conversion of the total angular and rate error to exponential 
form leads to the complex vector representation 

<w iAt -iAt , 

fixy = ^ + ifly = fixy + Bje + Bl+e (5 

where 


n t> _ -iAt 

fJjcy = Bj_e + B2© 



and the complex constants are given by 



- rr 4- -f I 4- 4. P_ 



c. Initial Error Contribution 


The total errors and, which correspond to the initial 

conditions ^q, 0y O) a, cp 0 , and 0 O will be considered first, A 
simple geometrical interpretation of these error traces is possible. 

For ftjjy this interpretation follows from the trace of the velocity- 
error in the XY body axis plane, as shown in figure 5. 


Y 



The path described by the tip of the vector is an ellipse in the 

body- fixed plane. 


The characteristics of this ellipse are derived from an examina- 
tion of (52) and (53). The semiaxes a and b are determined as 

r=r~ /(^x^xo + (^y^Vo /^y) 


^ = 1 y?y sj^x^x.0 7 ® 


(57) 
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and the angular position of the rate error vector is given by 



The quadrant for the angles corresponding to the inverse trigonometric 
functions tan -1 ( ) in ( 58 ) and in all subsequent equations is deter- 
mined by the sign of the numerator and denominator of the term in the 
brackets. When both numerator and denominator are positive, the angle 
is in the first quadrant; when the numerator is positive and the 
denominator is negative, the angle falls in the second quadrant; when 
both numerator and denominator are negative, the angle falls in the 
third quadrant; and when the numerator is negative and the denominator 
is positive, the angle falls in the fourth quadrant. 

The position of the major axis of the ellipse is determined by 
the relative magnitude of I x and Iy . If ly > T x> then the major 
axis coincides vith the Y body axis and the maximum angular rate occurs* 
about this axis. Conversely, if lx > Iy, then the major axis and the 
maximum angular rate lie along the X body axis. The period of revolu- 
tion is — and the vector rotates in the direction of the 

precession rate A. When I z is a maximum inertia, this rotation is 
in the direction of spin; when I z is a minimum inertia, the rotation 
is against the direction of spin. 

The trace of the rate error vector can be directly compared with 
the results of Poinsot’s geometric construction (ref. 12), in which the 
path of the instantaneous rate vector on the ellipsoid of inertia is 


called the polhode. For the present solution the rate vector is 
restrained to move in a plane normal to the Z axis, which is a principal 
axis of the inertia ellipsoid. The polhode projection onto this plane 
has been developed by Thomson (ref. 6, page 124) and yields a curve 
whose shape is defined by the relation 

^x(^oAc) + ^y(ly%-) = ^jc^x^xo) + ^y(ly^yo) (59) 

This relation describes an ellipse, with semiaxes given by (57). Since 
the polhode projection is proportional to the rate vector trace derived 
in this analysis, the approximate solution will exactly represent the 
spacecraft rates when the variation in the spin rate is negligible. 

The angular trace with respect to the Xj and Yj axes is illus- 
trated in figure 6. 



Figure 6.- Attitude error trace for initial conditions. 
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The path described by the tip of the Sj vector is generated by a point 

moving on a displaced ellipse, which in turn is rotating at the spin 

rate. From (55) and ( 56 ), the center of the moving ellipse is located 

by the vector sum of the initial attitude error a,-, and of the initial 

angular momentum ratio term i ~^ x — - - * ■. VPV° 1 . The radius a shown 

L oi z J 

on the figure is thus 




(ix^xo) + (^y^Vo) 


( 60 ) 


while the semiaxes of the rotating ellipse become 


b = 


CT -^z V 'Sc 


1 x^X 0 pT x f + (iy^yo v/^y) 


C = 


Olz J Ay 


. J( Ix^xo + (%%o >/^y) 


(61) 


and the precession of the attitude error vector within the ellipse is 
specified by the angle 




tan 


-1 



' I y%o 
,-^-x^xo J ^x > 


(62) 


When A is rational, the path of the attitude error trace is closed and 
has a period of 2kjt, where k is the least common denominator of a 
and A. 


The trace of the attitude error vector is in agreement with the 
general properties predicted by MacMillan (ref. 12) for the torque-free 
motion of a rigid body with respect to a unit reference sphere. This 
sphere was drawn about the fixed point of the spinning body as a center} 
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and the motion of the body Z axis about the fixed-momentum axis was then 
described by the trace of the Z axis on the unit sphere. The vector 
trace, introduced in the present analysis, can be considered as the pro- 
jection of this Z axis trace onto a plane perpendicular to the Zj axis. 

It should be apparent that the Zj axis which is arbitrarily 
defined as the fixed space axis corresponding to the initial position 
of the Z axis, will not generally coincide with the fixed-momentum axis. 
By assumption, however, the angle between these two axes is small. Hence, 
the shape of the traces about the fixed-momentum axis should be approxi- 
mately retained in the plane normal to the Zj axis. The fixed-momentum 
axis will appear as a displaced point on this attitude error plane. 

In figure 6, the fixed-momentum axis projects as the center of 
the rotating ellipse. The attitude error oscillates between two con- 
centric circles drawn about the ellipse center. The radii of these 
circles are given by the minor and major semiaxis of the ellipse. The 
similarity of this motion with that depicted in figure 6l of MacMillan's 
test is obvious. 

Upper bounds of the values for the rate and attitude error 
magnitudes, as developed from (45) and (46), are 
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and 


a 


lim 


- /(lx%o ) 2 + (^y^-o ) 2 + -~£- 


-1- + -i- 






J^x^xo /£> + (^y^yo 


( 64 ) 


The rate limit (63) gives the major semiaxis of the rate error trace and 
is equal to the max imum rate error. The attitude limit ( 64 ) corresponds 
to the sum of the center radii and the major semiaxis of the attitude 
ellipse, and will be greater than or equal to the maximum attitude 
error. 

Several interesting trends may be observed from the geometrical 
development and the relations for the upper limits of the errors. When 
the spacecraft inertia I x (or I y ) approaches I z while the second 
inertia I y (or I x ) remains different from I z , then the rate and 
attitude ellipses become very elongated. Small rate errors induced 
about the second inertia axis by impulsive torques or other disturbances 
can thus lead to large total attitude and rate errors. An example is a 
cylindrical configuration spinning about an axis normal to the axis of 
symmetry. 

When the spacecraft inertia I x (or I y ) is very much larger 
than I z , excessive attitude errors will be produced by s m all body rates 
and tumbling may occur. This result, however, is not surprising since 
the inplane angular momentum is now much larger than the spin momentum. 
Examples here are slender cylindrical satellites and missiles spinning 


about a minimum axis of inertia. 
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One may note that the smallest errors vill be produced when both 
I x and Iy are much smaller than I z , and the spacecraft configuration 
approaches that of a disk. 

The contributions of the errors and a to the limiting 

errors and |a| lim for a given disturbance will be omitted 

in the remainder of the analysis to avoid undue complications of the 

limiting error relations. These error terms could, however, be readily 

included if this is desirable for a particular disturbance. 

d. Impulsive Torque Contribution 

If the initial error terms are taken as zero, the total errors 

for the impulsive torques are equivalent to those for the initial rate 

error terms. The geometrical representation and the maximum error values 

for the vectors corresponding to the initial errors will thus hold for 

T 

the impulsive torques if n xo is replaced by 0y Q is replaced 

T I x J 

by and cp Q and 0 O are taken as zero in ( 57 ) - (64) and in 

figures 5 and 6. 

2. Step Torques 

a. Time Histories 

The spacecraft attitude control system and external sources, 
such as gravity gradients, may also exert torques about the body axes. 

For this example, consideration will be given to constant step torques 
of the form 


M s M x + iM y = (T x + iTy)u(t) 


(65) 
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and. the associated forcing function 


P --L 
lx 


T x 6(t) - AJT U(t) + T/(t) + \T x U(t) 

L J I yl_ 


The hody rate error is found from (2l), (U2), and (66) by sub- 
stitution of the solution functions of table 1, and is 

Sin “ cos ^0 

+ — T v sin At + -£~(l - cos At) / (£ 

T yL^ * JJ 


The attitude error a is determined in a like manner as 


c, = a + 




- cos At 1 - i cos at 


+ T v ^ sin At - — sin at 
y A a 


+ i S T y Twfe " cosAt ) -^ cosrt 


+ T y — sin At sin at 


The initial error contributions 


b. Total Errors 


and a are defined as before. 


A transformation of these components to polar form yields 


~ iAt -iAt 

fijjy — fijjy + Bj + B^e + B^e 


( 69 ) 
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(71) 


with 
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The trace of the total rate vector due only to the step torques is shown 
in figure 7* 



Figure 7 -- Rate error trace for step torques. 

This trace is now an offset ellipse which intersects the origin at time 
zero. 

The elements of the rate ellipse can be determined from ( 69 ) 
and (70) . The center of the ellipse is located by the radius 

a = 
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The semiaxes of the ellipse are found from 


b = 


3L. 


\fht V \ \ 


c=^ 


2 2 

T Z + v 


I y vPVv^V 




and the angular location of the rate error vector is 


$-u = tan~ \ £ tan 


At - tan' 


-l f T X J\ 

\ T y \PV- 


( 75 ) 


( 76 ) 


The position of the major semiaxis of the ellipse is again dependent on 
the relative magnitudes of I x and I y . If I y > I x , the maximu m semi- 
axis is parallel to the Y axisj if I y < I x , the maximum semiaxis is 
parallel to the X axis. The motion of the ft vector is in the direc- 
tion of A and has a period of 

An extension of Poinsot's development (ref. 12) to the motion of a 
rigid body under step torques appears possible. The fixed reference 
point, with respect to which the polhode projection is generated, lies 
along the maximum angular momentum vector possible for the body. The 
angular accelerations vanish for steady spin about this axis of maximum 
angular momentum. By referring to (2) and ( 3 ), one notes that the 
associated coordinates for the fixed point are proportional to 



for the constant step torques. 

The shape of the polhode projection corresponding to this fixed 

point is defined by (59) - Substitution of (77) into (59) yields the curve 

2 2 

yw 2 + vwi 2 = + ?- (78) 

An inspection of figure 7 shows that the polhode projection is indeed 
represented by the rate vector trace. The fixed point coincides with the 
center of the trace ellipse and the equation of the ellipse becomes ( 78 ). 

The inertial attitude error corresponding to the step torques 
yields the trace shown in figure 8. 



Figure 8.- Attitude error trace for step torques. 
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This trace results from a point moving along an ellipse, which remains 
fixed with respect to the rotating radius of a stationary displaced 
circle. The motion begins at the origin at time zero. 

The center of the stationary circle is located by the radial 
coordinate 


a= ^f^ + T - 


2 . m 2 
y 


( 79 ) 


corresponding to the ratio of the torque to the spin vis viva, 
radius of the stationary circle is 


The 


. i Kf + /if 

a VVW/ V'xXx.) 


and the center of the moving ellipse has the angular position 


♦b - fft + 

yx^-x^xj y T x j 


(8o) 


(81) 


The semiaxes of the ellipse are computed from the parameters 


c = 


i K 2 , v 

al z ** 


d = 


JT X 2 . ^2 


and the semiaxis 


aI z V V \ 

I c | makes the angle 


~\ 


(82) 


c 2 texj 


( 83 ) 
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with the radius h of the stationary circle. The angular coordinate of 
the tip of the attitude error vector within the ellipse 



completes the development of a point on the error trace curve. 

The motion is a closed curve when A and a are both rational. 
The corresponding period is given by 21a t, where k is the least common 
denominator of a and A. 

The general properties of the motion can be readily interpreted 
from figure 8, if one recalls that this figure represents the projection 
of the Z axis trace onto a plane normal to the Zj axis. The fixed space 
axis corresponding to the axis of maximum angular momentum projects as 
the center of the stationary circle. The motion of the body is bounded 
by two circles, concentric with the stationary circle and tangent to the 
moving ellipse in the figure. The outer circle represents the msy-ii-mmi 
Z axis excursion relative to the fixed momentum vector and can only be 
approached from the inside. The inner circle represents the minimum 
Z axis excursion and is approached from the outside when the ellipse does 
not contain the fixed momentum reference point. When the ellipse con- 
tains the momentum reference point, then the inner circle is crossed by 
the Z axis trace and is approached from the inside. 

Upper bounds of the error magnitudes are found from (45) and 
(46), with the result 
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and 




( 85 ) 



These upper bounds consist of the sums of the radial vector magnitudes 
and the semimajor axis of the error ellipse, and will obviously be 
greater than or equal to the maximum error values. 

For cylindrical spacecraft spinning about an axis normal to the 
cylinder axis, large rate and attitude errors will be produced by torques 
applied about the cylinder axis. Conversely, torques applied about the 
normal inertia axis in the spin plane will have little effect on the 
spacecraft motion. 

Near-cylindrical spacecraft spinning about a minimum inertia axis, 
so that I z is much less than I x (or Iy), will now be stable if the 

p 

applied torques do not approach the spin vis viva term I z a . 

5. Step Products of Inertia 
a. Time Histories 

Crew or cargo movements within the spacecraft may be represented 
by equivalent masses mj with variable position coordinates xj, y^. 
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and Zj. The movements of the equivalent masses fall into two categories. 
The first of these includes arbitrary nonperiodic motions along linear 
paths to some final position. From previous results for symmetrical 
spacecraft, it appears that the largest rate and attitude errors for 
such a motion are less than or equal to those for instantaneous motion 
to the final position. The introduction of step products of inertia 
corresponding to the final position coordinates of the moving masses 
gives a limiting case for this type of motion. 

The coordinates of the jth mass may thus be written as 


x j " x jo 
= y jo 

z j = ^ 0 U(t) > 


and the corresponding forcing function is 




T \ XZ 

x y ' 


5(t) + 0%(t) 
B(t) + a%(t) 


+ V-xz^W + a%(t)J 

- 7 Sc X yz jl.(t) + a 2 U( t)] 1 


( 87 ) 


( 88 ) 


where the products of inertia now take on the constant values 

n 


■'•xz ~ / m j x ,)o z jo " m s x so z so 
n 

•'■yz = 7 m j y j o z jo “ m s y so z so 




( 89 ) 
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The solution for the rate error becomes 


V =‘ s xy + (° 2 “ 


lx v/\ 


Ixz La At - -J? 
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2 -*J /V 


~r sin Xt 


x y /fy 


i^-La At - -HL_\ + *£_ sin At 


/x^\ (T 2 - X 2 / ^ 


. (is . HsWt) 


hr/ 


(90) 


and the attitude relation yields 


a = a + 1 


^ \ /l'V'7. _ III 


, — - p^,.- 7 === cos At + yL: sin At] 

^yvV IxvWlv^ v^x 


; ) + 1 ( 


(ix/*T 




X V^-yz . *xz , . , 

— sin Xt - — eos Xt I 


!y y ^y/V /s 


. 0 (is. . Hii) 
/ W*y ^x-’-x J 


(91) 


b. Total Errors 

The vectorial representation of the total errors reduces to 


„ iXt -iXt 

V = “xy + ®3 + B 4 e + B 5 e 


(92) 




-Al- 



and 

~ i(a+A)t i(o-A)t i<jt , , x 

a I = aj + C^e + C^e + Cge (9^) 


where 



The trace of the rate error vector is illustrated in figure 9. 
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Once more one obtains a displaced ellipse. However, the ellipse now 
does not intersect the origin. 

The radius to the center of this ellipse is 


a 




and the ellipse characteristics are found from 







J 


(96) 


( 97 ) 


and 


*b 


tan 



( 98 ) 


The major semiaxis of the trace ellipse is parallel to the X axis when 
I x > Iy and is parallel to the T axis when Iy > I x . However, the body- 
axis with the largest rate error is determined primarily by the location 
of the center of the trace ellipse. The period of the counterclockwise 
motion is again the precession period an< ^ the direction of motion 

is in the direction of the precession rate A. 

To correlate figure 9 with Poinsot's development, note that the 
coordinates of the fixed reference point for the polhode are propor- 


tional to 
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and 


o* - - 1x2 

XI 

A x A x 


« 2 t 

o • = - " , hi 

yo 


Vy 


(99) 


from (2) and (3). This shift of the reference point from the coordinate 
origin to the principal body axis indicates that the maximum possible 
angular momentum vector lies along the new extreme inertia axis. Steady 
spin, for which the polhode reduces to a point, is thus possible only 
about the new maximum or minimum principal axis of inertia. In figure 9, 
this principal axis passes through the center of the rate ellipse, as 
specified by ( 9 ^). 

The interpretation of the polhode projection about the fixed 
reference point becomes scanevhat more difficult. Two terms now con- 
tribute to the polhode, namely the rotation of the extreme principal 
axis and the effective acceleration torque produced by introduction of 
the step product of inertia. The rotation of the extreme principal axis 
yields the initial rates given in (99) • The step introduction of the 
product of inertia yields the additional rate terms 


and 


n « 

XO I * 


o .« _ 

V 


(ioo) 


The polhode projection (39) corresponds to the rate vector trace in 
figure 9 > if the reference body rates and ^- 0 include both ( 99 ) 

and (lOO). 
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The angular position trace for the step products of inertia is 
presented in figure 10, 


Figure 10.- Attitude error trace for step inertia products. 



The trace is produced by a point on an ellipse, which remains fixed with 
respect to the rotating radius of a stationary circle centered at the 
origin of the inertial axis system. 

The radius of the stationary circle is 


a 
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and the ellipse center is defined by the angle 


*a 
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The ellipse cemiaxes are found from 


b 




c 








J 


( 103 ) 


and the angle between the stationary circle radius a and the ellipse 
semiaxis | h | becomes 


-l/^xz^yly' 


i Iyz A X I X j 


= tan 

The position angle of the attitude error vector tip is now 


(10’f) 


$ c = tan’"'*' ^ ^ tan 


At - tan 



(105) 


with respect to the ellipse semiaxis |b|. 

The period of the precessional motion is 2k it, with k taken 
as the least common denominator of a and A for rational a and A. 

The general properties of motion for the step inertia product 
are similar to those for the step torques. However, the fixed space 
axis corresponding to the maximum angular momentum vector projects as 
the origin of the XjYj plane, and the direction of the total angular 
momentum vector for the spacecraft is not changed in inertial space 
during the step crew motion. The spacecraft motion is bounded by two 
circles, drawn with center at the origin and tangent to the moving 
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ellipse. The nature of these boundaries has been discussed in the step 
torque analysis. 

Upper bounds of the errors are 



for the attitude error vector. 

As before, cylindrical spacecraft spinning about a normal axis 
lead to large errors for small products of inertia in the plane cor- 
responding to the two large inertias. Seme differences in the response 
for the step inertia products and that for step torques can, however, 
arise after the removal of the disturbance. For the step torques this 
removal can occur when the rate error vector passes through the origin 
of the body axis system, so that the only residual error is a constant 
attitude error corresponding to the attitude at the time of disturbance 
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removal. This fact may be of use in the design of pure attitude control 
system for spacecraft -which use constant torque pulses to reorient the 
spacecraft. 

For step products of inertia, removal of the disturbance could 
also null the rate errors if done when the rate error vector passes 
through its initial position. In practice, the determination of this 
position does not appear feasible without very exact values for the 
spacecraft and disturbance characteristics. The elimination of the body 
rate errors by the timely removal of the product of inertia is a very 
complex task, and in general will lead to both residual rate and attitude 
errors . 

Variable Products of Inertia 

A second category of mass movements within the spacecraft involves 
periodic motions, such as mass transfer along a circumferential path. 

The uncontrolled motion of the spacecraft is now similar to that of a 
spring-mass system with a periodic forcing function. The amplitudes 
of the rate and attitude error are correspondingly multiplied by a 
magnification factor and resonance may occur for particular mass transfer 
rates. The determination of the effects of periodic mass motions within 
the spacecraft is thus an essential prerequisite to the analysis and 
selection of the spacecraft control system. 

There are, unfortunately, an infinite number of possible periodic 
crew motions. The best approach to a study of these motions may be the 
formulation of a general forcing function, which permits the development 
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of stability criteria and rate and attitude errors for a number of 
representative periodic motions. Particular time histories can then 
be developed for special cases of the general forcing function. 

To arrive at such a general forcing function, assume that the 
periodic mass motions involve transfer of a single mass and take place 
in a spin plane perpendicular to the Z axis or along a line parallel to 
the Z axis. If one adopts a sinusoidal variation of the associated 
mass position coordinates, the complex forcing function takes the form 

F = Eq cos pt + Ej sin pt + Ep&Ct) + E^6(t) + E^b’Ct) (108) 

The coefficients Ej are complex constants which must be determined 
for particular mass motions, and p is the period of the motions. 

A number of characteristic motions included in equation (108) are 
listed in table 4 along -with the corresponding coefficients Ej. As 
seen from the table, the forcing function equation (108) comprises linear 
periodic motions parallel to each spacecraft axis in the spacecraft 
reference planes XZ and YZ and circumferential motion in a plane per- 
pendicular to the spacecraft Z axis. Other mass motions can be con- 
structed by combining the forcing functions in the table and by adding 
the forcing function equation (88) for the static products of inertia 
•with appropriate values of I x2 and Iy Z to the result. Any linear 
oscillation in a spin plane perpendicular to the Z axis or along a line 
parallel to the Z axis can be developed by this method. 

Only the motions described in table k will be considered as examples 
for the present analysis j the results obtained for these motions can be 
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readily applied to more complex motions obtained by linear combinations 
of the forcing functions given by the table, 
a. Stability Trends 

An assessment of the stability trends for the variable products 
of inertia can be obtained by examining the forcing function terms in 
equation (108) together with the solution functions of tables 1 and 3 • 
The uncontrolled spacecraft motion will be unstable when the roots of 
the governing equations (29) and (3*0 contain real, positive, nonzero 
terms or when the forcing functions produce resonance conditions. From 
table 1, it is apparent that the solution functions for S(t), S(t), 

and 5(t) contain only constant and periodic terns. These terms thus 
cannot cause divergence of the rate and attitude errors. 

The solution functions for cos pt and sin pt can, however, 
contain divergent terms and may lead to continuously increasing errors 
for special frequencies of the periodic motions. These special fre- 
quencies are |p| = |a| and |p| = | A | . In the first case, the rate 
errors remain bounded for all finite values of Eq and Eq, but 
precession of the spacecraft may result unless 


when p = 0 and 


Eq + iEq = 0 


Eq - iEq = 0 


(109) 


when p = - 0 . Since equation (109) holds for all the values of Eq and 
Eq in table 4, this precession will not occur for the examples considered 


here. 
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In the second case, when |p| = |A|, both the rate and attitude 
errors will tend to diverge for nonzero values of Eq and Ej_. Mass 
motions with this period thus exhibit definite trends toward instability 
and should be avoided. 

There is, of course, one other instability that may occur for 
the present solution. Fran the governing equation (29) for the body 
rates, one notes that the rate error will diverge vhen A 2 < 0, so that 
I x > I z > Iy or ^-y ^ ^ J x' This condition results when the Z axis 

is an intermediate axis of inertia and agrees with the well-known fact 
that the undamped spacecraft spin is stable only if the spin axis is an 
axis of maximum or minimum inertia. In terns of the moving mass param- 
eters, one may thus write 
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U=i J 
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( 111 ) 


as the alternative conditions for instability corresponding to the 
assumptions of this analysis. 
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The conditions in this section will serve t.n indicate possible 
instabilities for the rotating spacecraft. Since the rate and attitude 
errors for these instabilities will rapidly exceed the small angle and 
rate assumptions, time histories for these motions will not be discussed. 
If the unstable motions do occur, exact computer solutions should be 
used to assess their effects, 
b. Time Histories 

From equation (108) and the solution functions of tables 1 and 3, 
the rate error becomes 


o * / E 0 \ /Ei 

^xy = ^xy + ( To 3 l cos P* + ( — ol sin pt 


VA 2 - p 2j 


^ - v 21 


+ E]^5(t) + E^ cos At + Eg sin At 


(112) 


where 


e 5 = e 3 - 


E, 
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U 2 - p 2 / 


E6 ■ x 


(E 2 - A%) - 2 


while the attitude error is 


a = a + 


U7 2 - p 2 / 




> (U5) 


fPBl " lo'Eo) A /PE 0 + iaEp) . 

r s pt * j sin pt 


\a 2 - A 2 ; 


(AEg - iaE^)cos At - (AE^ + icrEg)sin At 


+ E 7 (cos at - i sin at) 


(114) 
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where 


Ey - E^ - 


ABg - igErr pE 1 - irTE 0 


cr 2 - A 2 


(cr^ - P 2 )(A 2 - p 2 ) 


( 115 ) 


c. Total Errors 

The vectorial representation of these errors yields 

%y = %y + ^ + B 4 elAt + ^e" 1 ^ + B 6 e lpt + Bye" ipt (ll6) 


with 




= E]j5(t) 


b 4 » \ (=5 - m s ) 


B5 = \ (E 5 + IEg) 


l/E 0 - iEi> 

B 6 = p 


2 A 2 - p 2 


By^fl 


l/ E 0 + ^l' 


,A 2 - p 2 


J 


(HT) 


and, from (l8) 

otj = Oj + i 


+ C^e * ^ 


(a+A)t + C 6 e i(a “ A)t + C 7 e i(a+p)t + C ft e i(a - p)t 


Cge 


(118) 
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where 


C 4 = 


Bit Br< B ( : n 7 'N 


. V, jiv 

ii.]j + ___ + _ +. — + < 

ct + A tr-A cr + p 


P tT - p 


c 5 = - 


b 4 


a + A 


C 6 = - 




C 7 - 


c 8 


a - A 

b 6 

cj + p 
*1 

a - p 


(119) 


J 


The total error relations given by equations (ll2)-(ll9) may now be 
evaluated for particular motions by substitution of the corresponding 
coefficients Eq, Ej* and E 4 from table 4. Only the 

vectorial representations and the upper bounds of the errors will 
usually have to be considered in an assessment of the effects of mass 
motion on the spacecraft motion. 

To determine the vectorial traces one needs values for E 4 , B 4 , 

B6, and 3j. These egressions have been developed and are shown 
in tables 4 and 5 for the motions described in table 4. The polar plots 
of the rate errors now follow directly from equation (ll6) and are 
graphed in figure 11 for a number of typical motions. 
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All of these traces are generated by a point on an ellipse, whose center 
moves along a second ellipse centered at the origin of the body axis 
system. Figure 12 shows this development of the rate vector trace. 


Y 



Figure 12.- Rate error trace development for 
periodic inertia products. 


The direction of the ellipse semiaxes | a | and | c | is first located 
by the relation 



p 

a > 0 


a 2 < 0 


( 120 ) 


The position angles for the fixed and moving ellipses then become 


= tan' 


-1 


and 


$ c = tan 
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- tan pt 
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The fixed ellipse has semiaxes |a| and lb I, and the moving ellipse 
has semiaxes |c| and |d|. The constants a, b, c, and d are 
determined from (ll6) and (llT) as 


a = Bg + Bj 
b = B 6 - By 

> 

c = + 
d = B 4 - B^ ^ 


( 122 ) 


where the values of the Bj terms are taken from table 5 • 

The direction of motion is determined by the signs of p and 
A, and the rate vector describes a closed curve when p and A are 
both rational. The period of motion is 21ot, where k is the least 
common denominator of p and A. 

The rate error trace results frcm the oscillation of the mass 
with respect to the geometric body axes (p-ellipse), added to the pre- 
cession of the geometric body axes about the principal body axes 
(A-ellipse). The fixed reference point for the trace falls on the 
origin of the geometric body axis system since the mass oscillation 
takes place about the origin. 

Characteristic traces of the attitude errors for the periodic 
motions are illustrated in figure 13 - 


(a) Circumferential mass motion. o>) X-axis mass oscillation in 



(c) Z-axis mass oscillation in XZ plane. 

Figure 13.- Attitude error traces for periodic inertia products. 
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These traces are produced by two ellipses which rotate at the spin rate, 
as sketched in figure l4. 



Figure 14.- Attitude error trace development for 
periodic inertia products. 


The center of rotation is defined by the coordinates 

B 4 ^ A B 6 , B 1 


and 


a = -iEk + 
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The direction of the ellipse semiaxes |b| and |d| is located by the 
angle crt, measured from the positive Y-j- axis when $ a = 0 or n and 
measured from the negative Xj axis when $ a = ^ or The angular 

coordinates within the ellipses now become 


and 



(125) 


while the ellipse semimajor axes |b|, | c |, Id I, and |e| are derived 

from 


b » - 


S6 + ^ 


V* + P ° - V J 


'\ 


and 


c 



P 



/15l_ + 

\a + A o - A / 





J 


(12 6) 


The motion is a closed curve when a, p, and A are all rational. The 
precession period is 2fcrt where k is the least common denominator 
of c, p, and A. 

The attitude trace for the spacecraft with periodic inertia 
products exhibits a change in the direction of the angular momentum 
vector in inertial space. This rotation of the momentum vector to the 
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center of the p ellipse results from the initial mass acceleration 
terms, which exert a torque on the spacecraft. The motion, in body 
coordinates referred to the fixed momentum reference point, can be 
visualized as the sum of the mass oscillation with respect to the body 
axes (p-ellipse) and the precession of the body axes about the principal 
axes (A-ellipse). In the intermediate inertial coordinates, this motion 
is rotated through the angle ot. 

Values for the upper bounds of the rate and attitude errors may 
be calculated from 


l^xyllim ~ l Bi J + 1^1 + l^l + l^l 

and 


Bl^ Br Bg By 
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(128) 


by using the coefficient equations for E4 and Bj from tables 4 and 5. 

From an examination of these coefficient equations, one notes 
that the maximum error relations for small disturbance frequencies p 
yield the error limit terms corresponding to the introduction of step 
products of inertia. Both principal-axis-rotation and acceleration 
terms result for the circumferential motions; only principal-axis- 
rotation terms appear for the radial and vertical oscillations. A first 
estimate of the limiting errors for the periodic inertia products can 
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accordingly be obtained from the appropriate step inertia product terms, 
when the disturbance frequency | p | is much less than the precession 
frequency I ?s.| . 

As the disturbance frequency increases, the error limits also 
increase. As expected, divergence of the errors is predicted when 
I P | = |A|. For a further increase in the disturbance frequencies, the 
error limits continue to increase. 

When the disturbance frequency |p| is much greater than the 
spin rate cr, the error relations for the periodic inertia products 
become directly proportional to |p|. Doubling of the disturbance 
frequency will thus double the resultant error limits, and large errors 
can be introduced by small, rapidly oscillating masses which may occur 
in onboard motors, generators, or other equipment. 

Several additional trends are indicated by tables 4 and 5« For 
circumferential motion, maximum errors result when the sign of the 
angular velocity p of the motion coincides with the sign of the pre- 
cession rate A. Motion at a negative spin rate (p = -a) will eliminate 
all but the initial acceleration effects. Motion at a negative or posi- 
tive spin rate (p = ±cr) also nulls the errors caused by vertical mass 
oscillations, but does not significantly affect the errors caused by 
the radial mass oscillations. 

D. Characteristic Disturbances for Symmetric Spacecraft 
A large number of spinning spacecraft will be symmetric about their 
spin axis, so that I = I x = Iy. These spacecraft include rockets and 
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ballistic missiles which are spin- stabilized to maintain their flight 
path angle under initial body rates) unmanned satellites which are spin- 
stabilized to maintain a fixed inertial position for communication and 
observation purposes) and large manned space stations which provide an 
artificial gravity field for their crew. A reduction of the general 
solutions developed in the previous chapter to the special case of 
symmetric spacecraft would accordingly have many applications. 

Most of the resultant solutions have been previously obtained by 
various approximation and numerical integration techniques and are 
scattered through the literature (refs, b- 20). The results in this 
chapter thus make no claim to originality, but do accomplish two impor- 
tant objectives. 

The first objective involves the determination of the form of the 
geometric error traces and of the maximum error limits for the various 
disturbances. These important properties of the motion have been only 
p amiably treated in the literature, and tend to be obscured by tne 
component form in which past solutions are primarily expressed. The 
simple trace geometry, that results from the complex vector representa- 
tion of the present analysis should be of considerable value. 

The second objective is the comparison of the approximate solution 
with applicable previous results. This comparison will point out the 
principal differences between the present and past solutions and will 
summarize the trends of the motion with variations in the disturbance 
and spacecraft characteristics. 
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The discussions of the polhode projections and of the motion 
representation by means of the unit sphere will not be repeated here. 
If desired, these relations can be readily deduced from the analysis 
for the nonsymmetric case and the equations of this chapter. 

1. Impulsive Torques 

For impulsive torque disturbances of the form 


M = M x + IMy = 05(t) = (T x + iTy)b(t) 


the forcing function (48) becomes 


F = |[l(t) + iA&(t]0 


(129) 


The total rate and attitude error equations (^9) and ( 50 ) yield 


fljjy = + •— ( cos At + i sin At) 


(130) 


and 


a = a + -5— ("(sin At - i cos At) + (sin at + i cos at)] (l3l) 
al" *- ~ J 


Initial error contributions in (130) and (l3l) are 


^xy = %yo( cos At + i sin At) 


( 132 ) 


and 


a = a Q (cos crt - i sin trt) + i^-hp^Q(cos at - i sin at) 


- (cos At + i sin At)] 
from (4o) and (4l) . 


( 133 ) 
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The total error vectors for the impulsive torques reduce to 

®xy = ^xy + (f) elAt C 1 ^) 

and 

a I = + i (^j 1 - e 1 ( 0+A ^ (135) 

vhere the initial condition vectors are 

^xy = ^xyo e * ( 136 ) 

and 

3^=00 + i^aj(l - e l(o+A)t ) (137) 

The rate vector trace for symmetrical spacecraft with initial rate 
and attitude errors is shown in figure 15. 



Figure 15.- Rate error trace for initial conditions 
and symmetric spacecraft. 



result is 
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(138) 


and 


l a lmax = - V + cos $ ) 2 + ( 0 o + ^xo + a sin ^ (139) 


where 
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From ( 139 ) ) it can be seen that the spacecraft inertia ratio -f- 

lz 

does not affect the maximum rate error. The inertia ratio does, however, 
enter into the relation for the maximum attitude error vhich decreases 
with a decrease in the inertia ratio. Spacecraft, whose inertia ratio 
approaches that of a flat disk, will thus yield the minimum total atti- 
tude error for a given set of initial error values. 

The error relations developed here have been partially described by 
Leon (ref. 5) for the case of a spinning symmetric rocket, and the 
applicable present results agree with his conclusions. His work does 
not, however, develop the detailed trace representation or the maximum 
error relations. Thomson (ref. 6) includes the trace representations 
for the initial errors, but his geometrical interpretation is incorrect 
(see page 20l). The initial error vector a 0 and the initial body rate 
are not generally orthogonal, as depicted in Thomson’ s work. 
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The error traces and maximum error values for the impulsive torques 
can be again found by setting ^ and 9 0 = e 0 = 0 in 

equations ( 136 ) - (l4o) and figures 15 and 1 6 . 

2. Step Torques 

For constant step torques described by 


M h Mx + iMy = TO(t) = (T x + iTy)u(t) 


the forcing function ( 66 ) yields 


F = |[|(t) - iA(j(tj] 


(l4l) 


for the symmetrical spacecraft. 

From ( 67 ) and ( 68 ), the error time histories are 
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The total error vectors are 
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Curves for the rate error trace are given in figure 17. 


Y 



Figure 17* - Rate error trace for step torques and 
symmetric spacecraft. 


The rate error trace is a displaced circle whose center is located by 
the vector The radius of the circle is and the period of the 

I Prf I 
A I 

Curves for the attitude error trace are depicted in figure 18. 
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radius 


y 


Qi ip T C~ 

X ^ — • The center of this displaced circle is located 

oAl 


by the vector - 


T 


a d I, 


and the period of motion is 2kn, vhere k is 
the least common denominator of cr and A for rational o and A. 
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Figure l 8 .- Attitude error trace for step torques and 
symmetric spacecraft. 


The maximum values of the errors become 

Ivl 


I max 


2 hx 2 + Ty 2 


I z - Ik 


(146) 


and 




i z + 1 


a L . = 

1 llim 


Iz - 1 

+ 1 


V/T X 2 + Ty 2 


(147) 


The error vector relations for the step torques and symmetric space- 
craft agree with Thomson’s results (ref. 6 , pages 198 and 207). Maximum 
attitude error limits for this disturbance have also been developed by 
Suddath and include the residual errors after removal of the disturbance 
(ref. l4, page 8 ). His limits, which are smaller than or equal to the 
limits obtainable during the torque application, do not represent the 
worst case and are thus somewhat misleading. 
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Note that the error limits are inversely proportional to the inertia 
difference |l z - l| and predict minimum errors for a disklike config- 
uration. When applied to a sphere the limits predict divergence, 
confirming the statement that a sphere cannot he spin-stabilized. 

3. Step Products of Inertia 

When step products of inertia of the form 


■ 7 . 

5=1 


^XZ “ / ^^-jo 2 ,)© “ m s x so z so 


n 

V 

Iy2 - / m^yj 0 Zj Q - m s y so z so 

j=i 

occur in symmetrical spacecraft, one may introduce 


*rz “ ^xz + iJ yz 


(148) 


to get the forcing function 


F = - 


L rz 


5(t) + ff 2 U(t) 


- 1 


5(t) + 




from. (88). 

The error solutions (90) and (9l) give 


+ P? 


2 2 > 


(149) 


(a - ^_^ cog kt + i sin At) - + i&(t) 


A 


( 150 ) 


and 




(a - A)(cos At + i sin At) 


3 


(151) 


for the symmetric spacecraft. 



In vector form, one has 




m is^lL ) e 1At . £. + „( t) 


°T = °I 


jae iCTt - Co - X)e 1(o+?l)t 


The polar graph of the rate error is illustrated in figure 19. 


a 

a * V xz yz 


v /V xz yz 


Rate error trace— » I 

Figure 19.- Rate error trace for step inertia products 
and symmetric spacecraft. 


The error trace is a displaced circle with center determined by the 

q 2j / 2 2 \ I 

vector - rz . The radius of the circle is ( - ~ A j Jl xz 2 + i y 
and the precession period of the motion is I — I. 
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The graph of the attitude error for the step product of inertia is 
given in figure 20. 



Figure 20.- Attitude error trace for step inertia products 
and symmetric spacecraft. 


The error trace now is produced by a point on a rotating circle with 

moves around the circumference of 

The 


radius ^ ^ ) J^xz + ^yz which 

a circle with center at the origin and radius (£)/ 


2 2 

I ^ + T ^ 
- L XZ x yz • 


period of motion is 21ot, where k is the least common denominator 
of o and A, and the trace curve does not close unless o and A are 
rational. 

Maximum error values become 




max 


2d 2 - A 2 


for the rate vector and 


a 


lim ~ 


AI 


2d - A 
AI 


/ 


*xz + *yz 


(15M 




xz 


+ I 


•yz 


(155) 


for the attitude vector. 
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The errors predicted by ( 152 ) to (155) do not agree with previous 
results (refs. 4, 6) for product of inertia disturbances. This might 

be expected, however, since the previous analyses have neglected the 
energy associated with the introduction of the product of inertia. 

The limiting errors (15*0 and (155) diverge for a spherical con- 
figuration. Attitude error bounds for a long slender cylinder (A -» -a) 
are three times as large as the attitude error bounds for a flat disk 
(A -> cr) . Rate error bounds for these limiting configurations are, 
however, equal. 

Variable Products of Inertia 

Forcing functions, stability criteria, and solutions for variable 
periodic products of inertia and symmetrical spacecraft retain the 
form (108) to (ll9). The symmetric spacecraft does, however, permit a 
simpler combination of the inertia products and yields a better under- 
standing of the effects of the spacecraft parameters on the error traces. 

Consider a product of inertia of the form 

Trz = ^xz + ilyz 

= Qz 0 (x 0 + iy 0 ) (156) 

as produced by the motion of a single mass. Circumferential, radial, 
and vertical oscillations which yield this inertia product can now be 
examined. 

a. Circumferential Mass Motion 

The circumferential mass motion begins at the point 
x Q = r Q cos pt 0 , y 0 = r 0 sin pt Q , z 0 and continues at a constant 
angular rate p around the perimeter of the spacecraft. 
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The associated position coordinates may be written as 




where 


x = r 0 cos p(t + t 0 ) 
y = r Q sin p(t + t D ) 
z = z 0 u(t 0 ) 

r o = R + y Q 2 


a yo 


j 


The solution coefficients are taken from table 5* with the result 


b 4= 3 

B 5 = - 

b 6 = - 

By = 0 


L rz 


(a + p) g - (A 2 - p 2 ) + p 
A - p 2_ 


rzp 

~2I~ 


I r z( a + p)‘ 

(A - p)l 


(158) 


J 


and 


Ct = 


pi 


rz 


2a - A 


a 2 . A 2 


3 


Cc - ~ 


L rz 

I 


(a + p) 2 - (A 2 - p 2 ) 


-\ 


(a + A) (A - p) 2(a + A) 


/ ( 159 ) 


c 6 = 


Cv = 


Plrz 


2(g - A)l 

(a + p)l rz 
(A - p)l 







Figure 21,- Rate error trace for circumferential mass motion 
and symmetric spacecraft f 
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is produced by an ellipse, whose center travels along the circumference 
of a circle. 

The radius of this circle is 



(162) 


and the ellipse center is determined by the angle pt, measured from the 
axis corresponding to the radial inertia product. The ellipse semiaxis 
I b | also lies along the radial inertia product axis, and the precession 
angle of the rate vector tip -within the ellipse is 


•b 



(163) 


vnere 


and 


c = 


(g + p)‘ 
A - p 


(g * p r 

A - p 


- (A + p) 


J 


~\ 


I c + I * 
X xz -“-yz 


) (iSh) 


- A 


/ ^xz + -*-yz 


define the ellipse semiaxes. 

The attitude error trace in inertial space, shown in figure 22, is 
derived by a rotating ellipse whose center translates along the cir- 
cumference of a displaced circle. 



(167) 


and 
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Rg + p) 2 * (x 2 - p 2 ) pX 

/i ^ + I ^ 
\j x xz T x yz 

L x-p g-x 

( 0 + A)l 

f(g + p) 2 - (X 2 * p 2 ) pg 

lx 2 + I 2 

V J-xz + ^yz 

L x-p g-xj 

(g + X)l y 


(168) 


respectively. 

Limiting error values can be computed from the upper bound relations 


and 


Kl 


lim 35 


(g + p) 2 - (A 2 » p 2 ) + p 
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(° + P)' 
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£ + T ' 2 
xz x yz 


(169) 


l a 


lim 


p(2g - X) 


a 2 -r X 2 


(g + p) 2 - (X 2 - p") 
(a + X)(X - p) 


V 




a; 


P 


cr + p 

2(0 - X) 


A - p 




T ^ 4 . T ^ 

x xz + x yz 


( 170 ) 


These results Indicate that the rate and attitude errors for the cir- 
cumferential motion -will be larger than those for the static products 
of inertia ^ for which p vanishes. The limiting errors increase with 
an increase in p and tend to diverge as p approaches A. For values 
of p greater than A, the error continues to increase with an increase 


in p. 
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For very small disturbance frequencies (|p| « |AI), the upper error 
bounds are given by 


and 


| ^xyllim 



2 .2 
2cr - A 

Axz 2 + lyz 2 

A 

I 

2o - A| \) 

q 2 2 

X xz ^ x yz 

A 1 

I 


(171) 


(172) 


As might be expected, the maximum errors now resemble the error 
bounds for the step products of inertia. The resultant bounds are 
essentially independent of the disturbance frequency. 

For very large disturbance frequencies ( p » a), one has 


IV! 


lim 


~ ^"1 P I J^-xz + 1 


■yZ 


(173) 


and 


a 


lim 


p(5cr - 3A) 

yixz + lyz _ 

0 2 - A 2 

I 


(17^) 


so that both error bounds increase linearly with |p| for this case. 

The relations (l7l) to (17*0 can be used to determine error limits 
and their variation with the disturbance characteristics, when the 
absolute disturbance frequency is much less than the absolute precession 
rate or when the absolute disturbance frequency is much greater than 
the spin rate. 
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Radial Mass Oscillation 


■\ 


The radial mass motion begins at the point 0, 0, z Q and is 

defined by the coordinates 

x = x 0 sin pt 
y = y Q sin pt 
z = z 0 U(t) 

J 




(173) 


The solution coefficients, taken from table 5, now become 


_ iplrz 
b 4 - — - 

-d __ ^P^rz 

% -~S~ 


1 O 2 + 2 Qh + p 2 

2 ,2 2 

A - p 


B/r = 


= i(g + P) I r z 
' 6 2(A - p)l 

„ _ -i(g - p) 2 I rz 
^ ~ 2(A + p)l 
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(176) 


and 


C4 = 
C 5 = 


c 6 = 



ip(2<r - A)l rz 
(o2 - A2)I 


iP*rz 

1 a 2 + 2ctA + p 2 

( o’ + A)l 

2 ' A 2 -p 2 J 


iP^rz 
2(d - A)I 

i(g + p)l r z 
2(A - p)l 


c 8 


- p)l r2 
2(A + p)l 






( 177 ) 
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Rate and attitude vectors are 



as determined from substitution of (176) and (177) into (ll6) and (ll8). 
The rate error trace is shown in figure 23- 

Y Y 



Figure 23.- Rate error trace for radial mass oscillation 
and symmetric spacecraft. 
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This trace is generated by two ellipse envelopes, whose semiaxes | a | 
and | c| lie along a line normal to the radial product of inertia axis. 
The characteristics of these ellipse envelopes are developed from 



(180) 


and the precession angles within th$ ellipses become 


$ a = tan | ^ tan pt 


is 

is 




$ c = tan I — tan At 


(181) 


j 


The attitude error trace is also derived from two ellipses and 


takes the form illustrated in figure 2b 



Figure 24.- Attitude error trace for radial mass motion oscillation 

and symmetric spacecraft. 

The ellipse semiaxes | b | and | d | initially coincide with the radia] 
inertia-product axis and precess at the spin rate. The center of the 
inner ellipse is given by 

„ _ p(2cf - A)~| \/lx z 2 + x yz 2 f ^ Qr 

* i — (1& 

and the semiaxes for the two rotating ellipses become 

fp(ff + A)"j yi XZ + lyz 1 
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d = 


e = 


a + A; 


pr 2 + 2cA + p 2 cr 

I 2 + I 2 
V J-xz T - L yz 

L A 2 - p 2 a - A^ 

I 

I a 2 + 2crA + p 2 A 

\/^xz + -'■yz 

L A 2 - p 2 0 - A 

I 

J 


from (179) • The precession angles within the ellipses are 

= tan"^^- tan"^ptj 


and 


= tan' 


tan”" 1- At 


The upper bounds of the error magnitudes become 

IVl: 
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21 


These error limits again increase as | p | increases and diverge 

approaches j A | . 

When the disturbance frequency is very small (I Pi « I A I] 

may use 


1 - 

cr 2 

llim _ 

A 


\/lxz 2 + 


l yz 


(lti3) 

(184) 

(185) 


(186) 
as |p| 
y one 

(187) 


L. 


108 - 


and. 
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0 

11 *= * 




xz 


+ X 


y Z 


(188) 


as the approximate upper bounds for the rate and attitude errors. These 
upper bounds correspond to the maximum principal-axis rotation for the 
inertia product and do not contain the disturbance frequency. 

When the disturbance frequency is very large (|p| » a), one has 


■w 2|p| JxJ + I, 2 


and 


\ *!xy Him 

I 


Mlim = 

2 p(2d - A) 

It 2 4. t 2 
V X xz " r x yz 

O 2 - a2 

I 


(189) 


(190) 


The upper bounds (189) and (190) increase linearly with jpj. 

A comparison of (l7l) to (l7^) with (187) to (190) shows that 
the errors for the radial mass oscillation will be smaller than those 
predicted for the circumferential mass motion, when the disturbance 
(p, I rz ) and spacecraft (a. A, i) characteristics are equal, 
c. Vertical Mass Oscillation 

The vertical mass oscillation comprises a periodic mass motion, 
which starts at the point x 0 , y 0 , 0 and is described by 

x = x D 

y - y 0 

z = z Q sin pt 


'N 

> ( 191 ) 
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The solution coefficients for this motion are 


and 
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(192) 


( 193 ) 


2(A + P )I 

The error vectors, as developed from (ll6) and (ll8), yield 


^ ^ - 1( ° g ' - (X * p). 1 ** + (X - pja-ip* 

2(A^ - p^)l L J 


( 19 ^) 


no - 



The rate error trace for the vertical periodic motions is given 
in figure 25. 



Figure 25.- Rate error trace for vertical mass oscination 
and symmetric spacecraft. 


This trace is generated hy a point on a circle, whose center traces out 
an ellipse in the XY plane. At time zero, both the ellipse semiaxis 
| a | and the circle radius vector lie along an axis, which is perpendic- 
ular to the mass position radius. The enipse semiaxes are computed 


from 
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The trace curve is produced by a circle which travels along the p 
ellipse. Both circle and ellipse rotate at the spin rate a, and are 
initially alined with the radial product of inertia axis. The center 
of the ellipse is determined from 

P \J- 


a = 


I XZ + Iy Z 


v a + A/ I 

and the ellipse semiaxes are given by 
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c = 

The center of the precessing circle is at 

-lft 

= tan — tan pt 

and the radius of the circle becomes 


(200) 
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(201) 


d = 


P(cr 2 - p 2 ) 

/ 2 2 

V Uz + 1 yz 

(u + A) (A 2 - p 2 ) 

I 


( 202 ) 


This radius moves with the precession rate A within the circular 

envelope. 

Values for the upper limits of the rate and attitude error mag- 
nitudes may be calculated by using 
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( 203 ) 


- 113 - 


and 
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1 + 


2 2 
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F77 


+ i| 
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cr - p 
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I 2 + I 2 
xz x yz 


(204) 


The upper bounds of the errors increase as p increases. The instabil- 
ity trend at |p| = |A| is obvious frcm (203) and (204). 

If | p | is small (| p| « | A|), these equations yield 




11m 


and 


a 


llim 


0 2 

hxz 2 + Iyz 2 
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a| \1 

r i* 2 2 + 


(205) 


(206) 


so that both error bounds are approximately independent of | p| and 
contain only the terms associated with the principal- axis rotation. 
If | p j is large (|p| » cr), the upper bounds reduce to 

■ ■ ~ 2 M f 1 ’ 


a 


xyllim 


^xz c + ^yz 


(207) 


and 


-llim = ^ § 


P| 7 1 : 


2 + I 2 

xz + x yz 


(208) 


as a first approximation. Both these error bounds are directly propor- 
tional to | pj . 

From (I7l)-(l7h), (187)-(190), and (2Q5)-(208) it can be seen 
that the error bounds predicted for the vertical mass oscillation closely 
resemble those developed for the radial mass oscillation. The attitude 


errors for the periodic vertical motions will, however, be somewhat 
lower than the attitude errors for the periodic radial motions, when 
the disturbance and spacecraft characteristics are equal. 

In concluding this discussion of the variable products of inertia 
for symmetric spacecraft, one observes that the variation of the l imi ting 
errors with inertia distribution is similar to that described for the 
step products of inertia. Flat disk configurations yield the minimum 
errors, and spherical configurations are unstable. However, the magni- 
tudes of the rate and attitude errors are now considerably larger than 
the errors produced by the step products of inertia. 

E. Controlled Spacecraft Characteristics 
1. Governing Equations 

The motion of the controlled spacecraft can be defined by a method 
similar to that for the uncontrolled case. The torques produced by the 
control system are now particular functions of the measured vehicle 
angular position and rate. These torques can thus be considered as 
forcing functions applied to the uncontrolled vehicle equations. The 
solutions of the resultant differential equations yield the spacecraft's 
angular position and rate errors, as before. 

The analysis begins with the selection of a control torque command 
and the development of the corresponding equations of motion. As an 
example, a linear control torque g will be introduced as 
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s + ig y I y ~ j^lx^x + k 2x J 


dis + + • • . 


+ i 


ft 

Klyfiy + K2y / % dt + K^yCp + 


■■■] 


( 209 ) 


where the Kj k are the physical control gains that must he provided by 
the stabilization system. The error signals are amplified by these 
gains, and a control moment, whose value is equal to the sum of the 
amplified error terms given by (209), is then applied to the spacecraft. 

Particular nonlinear control torques, which lead to governing 
equations of the form discussed in the literature (refs. 21-24) could 
also be considered. The present application will, however, be restricted 
to a discussion of the linear control functions in (209). These linear 
control laws can be readily mechanized and allow a simple interpretation 
of the mechanics of motion for various types of sensor inputs (such as 
those derived directly from rate gyros, stable platforms, or Euler 
angle computers). 

In an analysis of the Spacecraft's stability, it is easier to 
deal with nondimensional control gains. One may define such gains by 
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(210) 
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so that 


at g x and gy become 


g x - A ^ix^x + ^ ^ dt + (Jk^ x 9 + . . . 


gy = A k ly fiy + k 2y f ~ ” ~ j ^ dt + ak^yCp + 


The equations of motion can now be written in the form 
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f y j \% + a (^yz + X xz a ) + 
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® s ( z s y s 


- y s z s) 
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Zj(ff 2 xj + 2ayj - xj) + zjXj 
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z s (o 2 x s + 2ay s - x s ) + z g x s 


(214) 


with the spin rate a taken as constant and positive. 

The Euler angle relations may be expressed as 

q> - crG = fijj. 

0 + 09 = fty. 

and the simultaneous solution of (212) and (215) will specify the 
spacecraft's motion. 

2. Control Requirements 

In practice, the spacecraft's rate and inertial position errors must 
he kept within specified deadbands which are determined by the spacecraft 
mission requirements. Control of the spacecraft to these accuracies may 
be provided by a reaction jet system and a momentum storage system, and 
the torques that must be produced by the system actuators can be readily 
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determined if the required control torques about the body axes are 
specified. With a reaction Jet system, these torques are generated by 
variable-mass-flow or pulse-modulated Jets, and for a momentum storage 
system composed of control moment gyros or reaction wheels the corre- 
sponding torque components along the gyro gimbal or reaction wheel axis 
are computed and applied. Concern will be given to the actual mechaniza 
tion of such systems later; for the time being only the body axis torque 
necessary to stabilize the spacecraft will be developed. 

Since the magnitude of the angular error in inertial space is equal 
to the magnitude of the body- ref erred angular error, the damping of a, 
and will assure the adequate stabilization of the spacecraft with 

respect to both body-fixed and inertial frames. The problem is thus 
reduced to the determination of control torques that will damp a and 
P~xy to zero or to small steady state values. 

3 . Control Law Formulation 

To investigate particular control laws, the corresponding control 
torque functions g x and g y must first be specified. The associated 
governing equations follow from (212) and (215). The stability regions 
for the governing equations can then be defined by making use of the 
conditions developed by Routh (ref. ll) and Hurwitz (ref. 25). If the 
selected control torques allow stable solutions for the Euler angles or 
body rates, time histories and complex error solutions can be found by 
the Laplace transform technique or by numerical integration of the 
linearized equations. If no stable solutions* are possible, the control 
torques can be rejected immediately. 
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To illustrate the applications of this technique, a number of example 
control laws will be analyzed, 
a. Pure Rate Control Law 

For a pure rate control law, consider 


K 

g x = ^x = ^Ox^x 

x x 

K lv 

gy - ^y - ^Xy*V 






so that the moment equations yield 


and 


fi x “ ^Ix^ + 





(216) 


(217) 


fly - A^ly^y - 


¥*K = 


The Laplace formulation of the corresponding complex rate error is 

s 2 - ACk-^ + k-j^s + A 2 (l + k^k^jj^xyCs) = V(s) (2l8) 

where the transform of the effective forcing function is given by 



( 219 ) 
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The time solution of (2l8) consists of the sum of a general 
solution, for which f x (s) and f y (s) are set equal to zero in ( 219 ), 
and a particular solution of the complete equation (2l8), for which 
f x (s) and fy(s) are specified for the applied disturbances. The 
functions f x (t) and f y (t) are the explicit, continuous functions of 
time defined by ( 213 ) and (2ll), and do not contain the rate or attitude 
errors. 

The particular solution of (218) is directly dependent on the 
disturbance under consideration, but characteristic trends for this 
solution can be indicated when the general solution is a damped vibration. 
Step functions in (213) and (214) will lead to constant residual rate and 
attitude errors; impulse functions and their derivatives will lead to 
damped transient rate and attitude errors which approach zero as time 
increases; and sinusoidal forcing functions will produce residual sinu- 
soidal rate and attitude errors. The amplitude of these residual errors 
is reduced with an increase in the damping ratio. 

The actual development of time histories for the various applied 
disturbances will not be attempted here. The primary tasks of the control 
system are the minimization of errors during a disturbance and the 
elimination of residual errors after removal of the disturbance. Both 
these tasks can be accomplished by the selection of stable gains that 
yield large damping ratios consistent with realistic control systems, and 
do not require the development of time solutions. 
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Stable gains are gains for which all roots in the general solu- 
tion have negative real parts. The characteristic equation for this 
general solution is 


O 

s" - + k ly )s + A (l + = 0 


( 220 ) 


and stability of the complex rate error requires that 

•\ 


and 


-Ak iy > Ak^ 


k lx k ly > ’ 1 
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( 221 ) 


One should note that (220) has the form 


s + 2rj-p>jjS + — 0 


( 222 ) 


The damped natural frequency o>jj may be expressed as 


“to 


= I A | ,/ 1 + k-^k^ 


(223) 


while the damping ratio rp and time constant t D are 
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^ + k lx k ly 
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t D H 1 
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r rPN- Mk^ + k Iy) 


(225) 


respectively. 
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a desired 
yields the 


(226) 


(227) 


(228) 

for the specified real values of k-^ and k^. The stability condi- 
tions (221) are automatically satisfied by (226) and (227) for positive 
real values of the damping ratio and time constant. Control gains for 
particular damping characteristics can thus be determined directly from 
the above equations. 

The resulting regions of stability for the control gain functions 
k lx s 6 n A and k ly sgn A are illustrated in figure 27(a). The gain 
functions must be in the area bounded above by the rectangular hyperbolas 
k lx k ly = -1 and "the straight line k^ = - k ly . Of particular interest 
is the fact that either of the two gains can be zero. This means that 
damping of the spacecraft's rates is possible with torques applied about 
a single spacecraft axis and derived from a single rate gyro for that 


Selection of the control gains is generally based on 
time constant and damping ratio. For the present case, this 
relations 
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and 
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(a) Rate gains. 
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(b) Rate damping. 




(d) Rate integral gains for small (e) Attitude gains for small rate 
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Figure 27.- Stability characteristics for the controlled spacecraft 
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axis. Hence, a control system providing rate control torques about two 
axes has inherent redundancy in case of failure of one of the system 
torquers. 

An indication of the damping characteristics obtained for the 
pure rate control is given in figure 27(b). This figure shows the over- 
damped, critically damped, and damped vibration regions corresponding 
to the stable control gains. Critical damping occurs when r D = 1, and 
yields the conditions 


and 


k lx “ k ly + 2 


k lx = k ly " 2 J 


(229) 


from (224). Stable gains that fall outside of the straight lines 
defined by these conditions will yield overdamped spacecraft motion; 
gains that fall between the straight lines will yield vibrational, 
damped spacecraft motion. 

From (224) and (225), it is apparent that the time constant is 
smallest when the two rate gains have the same sign. To optimize both 
the time constant and the damping ratio, one of the two gains can be 
selected as zero. Single-axis rate control for the spacecraft should 
thus be quite efficient. 

The transformed complex position equation is 

(s + ic) s 2 - ^(h^x + k^y)s + A 2 (l + k -| y kiy^j^s) 

= V(s) + [s 2 - A(k| y + k}y) + A 2 (l + h^k^y) 


cm. 


( 250 ) 
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Since this equation contains a purely imaginary root, the solutions for 
the complex position error will be neutrally stable. The pure rate 
command is thus limited in its usefulness to those applications -where 
only rate damping is needed. An example of such an application would 
be a manned space station, where the functions of rate damping and 
attitude control are often provided by different subsystems and where 
the rate control law is used to command onboard momentum storage system. 


From (2l8), the rate error may be written as 



(232) 
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The response of the spacecraft Is analogous to that of a spring-mass 
system with forced vibration and damping. The free vibration term 
corresponding to the initial conditions approaches zero as time increases, 
and the forced vibration term corresponding to the applied torque and 
mass motion effects is multiplied by a magnification factor or is damped 
to zero. The magnification factor is a function of the control gains 
and decreases as the terms and kqy take on larger stable values. 

The attitude error for the pure rate control becomes 



where s^ and S 2 are given by (232). 

From inspection of (233), it is apparent that the initial errors Oq 
and ^xyo contribute both free and damped vibration terms to the 
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spacecraft’ s motion. Since the damped terms will vanish for large time, 
only the purely oscillatory terms need be considered. The amplitude of 
the motion due to the initial position error a 0 is not affected by the 
pure rate control law, as one would expect from the results of the sta- 
bility analysis. The amplitude of the motion resulting from initial 
rate errors can, however, be reduced by proper choice of the control 
gains for the particular spacecraft under investigation. 

In summary, one may thus conclude that the pure rate control law 
is adequate for the damping of the spacecraft’s angular velocities. If 
no initial attitude errors or reorientation requirements exist for the 
spacecraft, it should also be possible to select control gains which 
will hold the spacecraft to small oscillations about its initial position 
in the presence of crew motions and other internal disturbances. The 
latter function is particularly important in spacecraft with solar cell 
panels or similar equipment, which must be approximately maintained in 
a given inertial direction. The effectiveness of the attitude hold mode 
for the rate control law should, however, be checked by substitution of 
"worst case" forcing functions for the spacecraft into (233). 
b. Rate Plus Rate Integral Control Law 

When the control torque is derived from both rate gyros and 
integrating rate gyros, one may take 


Sx 


= + K2x 



k 3jAc + k 2x| 


'W\ [ t 


Ay dt 


>(23^) 


. / ft \ 

A i \ ft 1 

+ a* at)-* 

k 2x*V + k 2yy^ j J ^ dt 
u J 


and 
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with the governing equations 


“ ^lx^ + " ^2x f dt J = 


*\ 


and 


/ (235) 


S - ^ly^y - m 


+ Akpy f ' tyr dl^ — f^ 


J 


The characteristic equation for the general solution for the rate error 
is given by 


s - A(k^x + ^ly) s ^ + ^^(l "*■ kqjjk^yJs^ + A^(kpy — kgj^Js — A^k2xk2y — 


(236) 


and the stability conditions are 


and 


-Xk ly > 
k lx k ly > _1 
~ k 2x k 2y > 0 


^27 > ^2x 




) (237) 


( k l* + k ly) k 2x k 2y > ( k 2y “ 

+ (k-^ 


k 2x^[( k 2y - k 2x) 

+ k ly )(l + k^ ly )j 


Here, the nondimensional gains kj^ are defined as in (212). 

The governing equation for the general solution for the position 
error yields the relation 
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s + la 


s - + k^y)s^ + V^(l + kjjjk^y) 8 ^ 


+ ^ 5 (^2y " k 2x^ s * ^ k 2x k 2y] a ( s ^ = 


0 


(258) 


Once again, the attitude equation contains a purely imaginary root, 
leading to neutral stability of the position error. The trends of the 
spacecrafts motion for the rate plus rate integral control lav are 
thus similar to those for the pure rate control lav. 

The stability regime for this control lav is shovn graphically 
in figures 27(a) and 27(c). The nondimens ional gains kjjc* ^ly> ^2x> 
and k2y must nov be selected to satisfy the conditions (237)* The 
first tvo of these conditions are identical to those for the pure rate 
control lav and are given by figure 27(a). The next tvo conditions lead 
to stable motion in the second quadrant of figure 27 (c), subject to the 
last restriction which represents a compatibility relation betveen the 
rate and rate integral gains. This compatibility condition yields a 
hyperbola with the equation 


k 

+ 


2y 


2 + (k^ + k ly r kgxk + k^" 


( k 2y “ k 2x) 


( k lx + k !y)(l + k lx k lyjj = 


0 


(239) 


as sketched in figure 27(c). The stable region in this figure is then 
the area between the upper segnent of the hyperbola and the k^ sgn A 
and k^ sgn A axes. 
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If the nondime n s i ona 1 rate gains are small in comparison with 
unity, as would he the case for most practical control systems, then 
(239) can be approximated by 

(kix + k ly ) = (k2x - k2y) (2ho) 

and the resultant stability characteristics are given in figure 27(d), 
The spacecraft’s motion is stable if the rate integral gains are 
selected from a triangular area of the fourth quadrant for a set of 
stable rate gains. 

By examining (235) and (236), one notes that the modified char- 
acteristic equation for the rate error reduces to a cubic equation when 
either of the two rate integral gains vanishes. This special case would 
occur during single-axis control of the spacecraft and is thus of 
particular interest. 

Stability restrictions are given by the standard rate gain 
restrictions and the relation 


0 > sgn A > (kj^ + k ly )(l.-f k lx k ly )sgn A 


(24l) 


The resultant gains will fall on the boundary of the stable region in 
figures 27(c) and 27 (d). Accordingly, stable control of the spacecraft 


is possible with single-axis rate plus rate integral commands. 
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Returning to the governing equations (235), one observes that 
numerical integration will be necessary to determine the rate time 
histories for particular values of the control gains. Some general 
conclusions can, however, be drawn for the case where the nondime ns ional 
rate gains are chosen to be considerably larger than the nondimensional 
rate integral gains. Inspection of (23 6) shows that the last two terms 
in this equation are now small during an initial transient period after 
a disturbance, when the rates are large and their integrals are small. 
The rate equation is thus approximately equal to (217) and the rate 
gains can be selected from (226) and (227) to yield the desired damping 
characteristics during the transient period. As steady state conditions 
are approached, the rate integral terms will predominate and the result- 
ant control torques will tend to eliminate any residual rate errors. 

The net effect of this law will be a reduction in the gain magnitudes 
since the high rate error gains will no longer be needed to reduce 
standoff errors in the steady state condition. Damping and attitude 
hold characteristics for the rate plus rate integral law should thus be 
quite efficient. 

c. Rate Plus Attitude Control Law 

Next, consider control torques developed from both rate and 
attitude errors and given by 


g x = ^-(Kqjcfix + K3 X 0) = A(k ljc D x + ak 5x 9) 

l x 


( 2 ^ 2 ) 


% - ^ K iy a y + K Jy' p) - * ‘’ k 3y' I> > 
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The resultant equations of motion are 


% - NyA * m 




£ly - a?vkj x 0 = f x 


and 


^x I x 1 


Qy — ™ I ■ I “ aAK^y^P — fy 


( 243 ) 


J 


leading to the general characteristic equation for the complex rate 


error 


s - A(k]_x + 


aA ( k Ix + k ly) 


k ly) + ^ "** + k lx k ly)j 

3 


s^ + a 


^ k 3y 2> hc I x k 3x%r I y^ 


k 


s + ct 2 A 2 | k^kjy + k^ - k ly | = 0 


(244) 


The associated Hurvitz stability criteria yield the conditions 


-Ahj_y > Ak^ 




k Ix k ly > 


1 + 


fefl 


2 2 

f k 3y k x-*-x k 3x ^y^y) 


> aA(kj_x + k ly) 


} (245) 


and 


f k 3y A X I X kjxAy!y\ A / k lx + k ly) 


> ;l 2 


0 - A ( 1 + k lx k ly ) 


.2 1 2 / 
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- A 2 (l + ki^ly^l + ^^^(l + k iy “ k lx^ 


for stability of the rate errors. 
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After substitution of (215) into (242), the characteristic equa- 
tion associated with the general solution for the attitude error becomes 


A 


[• 




- o 2 ~h 

n -* : ’ 

( k lx + k ly) 

+ a 2 * 2 

^ — 
1 + (k]^ - k^ x )(k^_y + ^3y) 

= 0 


( 246 ) 


Corresponding stability conditions are 




^lx^ly ^ 



_ SaA > (fix* *1 y , \ 

W lx / \I X + % - Iz/ 
( k Jjc " k 3x^ k ly + k 3y^ > _1 


> ( 247 ) 



for the attitude errors. 
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Both (p^p) and (pVf) must hold for stable error histories. In 
most cases the nondimensional rate and attitude gains are small in com- 
parison with unity, and the stability conditions given by these equations 
can be approximately represented by 

A 

> Akpx 


k lx k ly > - 1 


and 


0 > 


c 3x 


c 3y 


V k lx + k ly/bz / \ k lx + k-ly/bs 


> 1 


(248) 


The rate gain restrictions are now identical to those for the pure rate 
law and are thus shown in figure 27(a). The attitude gain restrictions, 
illustrated in figure 27(e), yield a stable region falling between the 
lines 


and 








c 3x 


^lx 


+ k 


iy> 


<lx 


+ k ly/\ I z / 


c 3y 


\ k lx + k ly) 


L JL\ - 


c 3x 


V k lx + k ly / -*-2 


3 - i 




( 249 ) 


It is worthwhile to note that either of the two attitude gains may be 
zero, and attitude and rate damping of the spacecraft is thus possible 
with torques applied about a single spacecraft axis. The error signals 
needed in the calculation of these required torques can be taken from 
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:i n ingle rut.e gyro and a oingle sun sensor mounted on that axis^ and 
the mechanization of such a control system appears to he very simple - 
As seen from and (2^6) ? the determination of the damped, 

rate and attitude errors again requires the numerical integration of (243). 
One observes, however, that both damping and reorientation control can he 
provided by the rate plus attitude control law, and that this law can 
correspondingly he used to maneuver the spacecraft. In comparison, the 
rate and the rate plus rate integral control laws were restricted to 
holding an already established inertial position. 

In addition, it may be noted that pure attitude control, for 
which k-^ = k ly = 0, will result in several zero coefficients in (244) 
and (246). The associated spacecraft motion is, at best, neutrally 
damped and may diverge for certain forcing functions. Pure attitude 
control then provides no damping of the spacecraft rate and attitude 
errors. This conclusion is, of course, in agreement with previous 
re suit s . 

d. General Considerations 

Other control laws may be investigated in an identical manner by 
selecting the control torques, developing the complex governing equations, 
and defining the resulting stability regions for the spacecraft's motion. 
If the control system gains are chosen to satisfy these criteria, the 
motion of the spacecraft will he damped. The determination of time 
histories of the controlled spacecraft motion becomes a rather tedious 
task, however, and is perhaps easiest if the governing equations (212) 
and (215) are programed on a digital or analog computer. Having assured 
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that thn selected control gains lead to stable motion of the space- 
craft, cases with particular disturbance functions can then be run on 
the computer to determine the spacecraft time histories. 

Extension of the method to include nonlinear control commands 
is possible, but nonlinear techniques (refs. 21-27) must then be used 
to define the stability of the governing equations. A preliminary 
selection of the type of on-off control commands, as represented by 
step torques, may be made by noting that the spacecraft response for 
an amplitude -limited control system with high gains approaches that 
for an on-off control system. Sensor inputs and the signs of the 
control torques may thus be chosen from the proportional analysis. 

Time histories for the on-off commands can then be obtained by substitu- 
tion of the corresponding step functions in (42). The solutions for 
the body rates and Euler angles now are found by a piecemeal process, 
and the forcing and solutions functions change whenever the deadbands 
for the on-off system are crossed ? 

Linear control optimization (ref. 28), as represented by a 
minimum mean square error criterion, may also be considered. Maximum 
torque or error limitations can be included as restraints in such an 
analysis. 
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h . Co n t^raL nystem Selection 

Having determined a control law which leads to acceptable damped 
motion of the spacecraft, one must next select physical systems that can 
develop the actuator-torque histories required by the control law. The 
choice of such control systems is generally made on a minimum- launch- 
weight basis; and relations between the control system weight and its 
impulse or momentum storage capacity are needed to evaluate the compara- 
tive merits of various control hardware. Preferably, these relations 
should not necessitate the detailed design and optimization of competitive 
systems for a particular spacecraft. 

An empirical representation of the total control- system launch-weight 
in terms of the angular momentum or impulse provided by the system will 
be used in this analysis. Such a representation gives reasonable 
approximate values for the control system launch weight, and allows the 
rapid comparison of different control actuation schemes. Furthermore, 
the empirical results are completely independent of the spacecraft 
inertia characteristics or dynamics. 

There remain then two tasks, namely sizing and implementation of 
the control system. To size the system, one first determines the space- 
crafts angular momentum envelope by integration of the torques cor- 
responding to simultaneous application of all "worst-case” disturbances. 
The launch weight for the control system can then be developed fran the 
empirical data, and a preferred system concept can be selected. The 
implement? it ion of this concept requires the solution of the control 


- 138 - 


system equations to define the actual torques the system must generate 
in order to provide the desired control law and eliminate cross coupling 
moments caused by any angular momentum stored in the control system. 

Control system components may he divided into the general classes of 
momentum storage units and reaction control units. Momentum storage units 
comprise reaction wheels, single-gimbaled control moment gyroscopes, and 
double-gimbaled control moment gyroscopes. Reaction control units consist 
of reaction jets with variable mass flow or pulse modulation. 

The angular momentum envelopes for these two classes of control 
components are then given by 


%1x + ^My 
-t d 


= I x fx + ifly 

- o J L u 


r d 


f y dt 


( 250 ) 


for the momentum storage units, and 


^ = %x + iH Ry 

- [u J** |fx|it] + ijly |%|dt 


( 251 ) 


for the reaction control units. The integration is carried on over t^, 
the time interval of application of the "worst-case" disturbances, 
a. Reaction Wheels 

For sizing purposes the reaction wheel will he taken as a flywheel 
which is accelerated by means of a torque motor to produce reaction 
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torques on the spacecraft. A sketch of a control system using two such 
wheels is shown in figure 28. 


z 



Figure 28.- Spacecraft control with reaction wheels. 


From manufacturer's data*, the basic weight of a reaction wheel 
(ref. 29) with a minimum alternating-current motor configuration is 


. Htj 

w w = 6.3 + 170 — 


(252) 


where Wy is the total reaction wheel and motor weight in pounds, Hy 
is the angular momentum capacity in ft lb sec, and u>g is the synchronous 
wheel speed in rad/sec. 

Similarly, the reaction wheel power may be derived from empirical 
data (ref. 30) as 


P W = 2 -T7 G^Og 


( 253 ) 


JL 


♦Reproduced and used by courtesy of the Bendix Corporation. 
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where Gg is the stall torque in ft lb. If one introduces a power 
weight conversion factor of a lb / watt , then the equivalent power system 
weight for the reaction wheel is 


Wp = 2.77 aGgcng 


(254) 


and the total weight chargeable to one reaction wheel becomes 


Wyip = Wy + Wp =6.3 + 170 + 2.77 aGgCDg 


To optimize the total weight for a given angular momentum and stall 
torque, one differentiates the total weight with respect to wheel speed 
and equates the result to zero. Substitution of the corresponding 
wheel speed into (255) gives 


W WT 


6.5 + 21 



(256) 


The control torques may be assumed to be sinusoidal with ampli- 
tude Gg and frequency A, so that one can take 


g = -G s e 


( 257 ) 


as a good approximation to the control moments. Since these control 
moments are equal to the total rate of change of the angular momentum 
components for the reaction wheels, it follows that 


Gg cos At = Hyjp — aHyg 
Gg sin At = H^2 


( 258 ) 
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or introducing 

H xy - %1 + 

(259) 

one obtains 

Crge*^ = H X y + ioHjjy 

(260) 

This expression 

can be integrated to give 



“ 

(261) 

for no initial wheel momentum. Maximization of (261) 

further yields the 

value 

11 _ °s 
% ' a + \ 

(262) 


for each of the two reaction wheels. 

Substitution of (262) into (256) leads to the expression 

= 6 . 


Wi 


WT 


>*•» 

>•3 + 21.7 \/a(cr 


+ A) 


% 


(263) 


which is plotted in figure 29. 



Figure 29-- Variation of total reaction wheel weight with 
required angular momentum. 
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The weight of the reaction wheel control system is now given by 

W T = + w WT(%y) (2 64) 

where Wyij(Hj^ x ) and Wwr(^My) are taken from the figure, using the 
values of and %y previously determined from (250). 

One should note that the spin rate a must be very small if 
reaction wheels are to be efficient. As an example, a power conversion 
factor of 1 lb/watt and a spin rate of 0.25 rad/sec for a flat disk con- 
figuration would yield Wrp = 3>07^ lb for an angular momentum requirement 
of 100 ft lb sec along each axis. Since such exorbitant weight penalties 
are impract ical, reaction wheels generally are inacceptable for the 
damping control of spinning spacecraft. 

If such wheels are used for spacecraft with very low spin rates, 
the governing torque relations become 

g = -G = ~Hxy - icrHxy (265) 

and 

T z = GyH x - GxH y (266) 

where G is the complex torque applied to the reaction wheels, and T z 
is the cross coupling torque applied to the spacecraft by the control 
system. Since the body rates are small, this sinusoidal cross coupling 
torque is relatively small and its effect on the spin rate a will be 
neglected. 
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The desired control torques are then obtained by directly applying 
the reaction wheel torques 


G - G W1 + iG^2 = -g (267) 

by means of the wheel actuators, 
b. Control Moment Gyroscopes 

A control moment gyroscope consists of a flywheel which spins at 
a constant speed and is mounted on a single or double gimbal arrangement. 
Control torques are now developed by precessing the flywheel. Torque 
actuators mounted on the gimbals provide the necessary precession torques. 
Sketches of control systems using single- and double- gimbaled gyros are 
given in figures 30(a) and 30(b), respectively. 
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(a) Single- gimbaled gyros. (b) Double-gimbaled gyro. 


Figure 30.- Spacecraft control with control moment gyros. 
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Since the weights of single- and double-gimbaled gyros do not 
differ appreciably, launch weights for both these units will be assumed 
to be identical. The basic weight of a control moment gyroscope can 
again be developed from manufacturer's data (ref. 31) and becomes 

W G = 1.37 Hq 0,68 (268) 

Power requirements now are derived primarily from the windage and 
friction. losses for the flywheel, and can be approximated by empirical 
data derived from computer analyses (ref. 32) as 

W p = 1.47 a %°* 362 (269) 

The power required by the gimbal actuators is small and will be neglected. 
The total weight of the gyro is then 

W GT = W G + Wp = 1.37 I^ 0,68 + 1.47 aHc 0 * 562 (270) 

Launch weights of the gyro are plotted against angular momentum in 
figure 31, and a comparison with figure 29 shows that the total weight 
for a system using gyros is much less than that for a system using 


reaction wheels. 
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Figure 31- - Variation of total gyro weight with, required 

angular momentum. 


In comparing single- and double-gimbaled control moment gyros, 
one observes that the weight for the single- gimbaled gyro system is 
derived from 


W T - + w GT( H My) 


(271) 


while the weight for the double-gimbaled gyro system becomes 


W T " W GT ! 


V %x 2 + %y 2 


cos a n 


(272) 


■where a,g is the limiting gimbal angle and and H^. are again 

found from (250). Since ag is generally 6o°, the double-gimbaled gyro 
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system is somewhat lighter than a control system using two single- 
gimbaled gyros. 

To derive the gimbal actuator commands for the implementation of 
the desired control laws, one notes that these torques are again derived 
from the total rate of change of the angular momentum vector for the 
gyro system. Thus, the complex control torque is 


g 





and the cross coupling moment becomes 


(273) 


T z - -H z + — G^Hy 


(274) 


The minor changes in the spin rate a due to T z will be neglected in 
this linear formulation, and the necessary control commands are now 
found by expressing E^y and H z in terms of the gimbal angles. 

For the single-gimbaled gyros this gives 


H xy = H x + iH y = Hqi sin 0 g - iH(j2 sin cp 






and 


> (275) 


H z = cos 0g + H G 2 cos cp g 


and 


The governing torque equations now reduce to 
g = -Gxy = - ||H G1 (e g + fiy)cos 0 g + Hq 2 (ct sin cp g + Gy cos cp g )j 
+ i |%2(^g “ P g + %i( 0 sin e g - cos 9 g )j^ 


(276) 


T z = ^(Og + fiy)sin 0 + HqpfPg - Gbc)sin cp 


(277) 
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where the small girnbal accelerations haye been neglected. Gimbal torques 
may be commanded directly from 

Gxy = G x + iG y = -8 (278) 

in an open- loop system. 

For the double-gimbaled gyro the angular momentum components 
along the spacecraft’s axes become 


h* = Hq( sin 9 - i sin q> g cos 9 g ) 


( 279 ) 


and 

H z = Hq cos 9g cos 9g 
The torque relations thus are 
g = — G xy = “( G x iGy cos cp g ) 

9 g + c sin cp g + Qy cos cpgjcos 9 g - i|(q>g + G x )cos 9 g cos cp g 




- (c + 0 g sin cp g )sin 9 g 


and 


3 } 


T z = -Gy sin cp g = % j( 9 g + G^sin cp g cos 9^ 
+ (^ + 9 g cos <p g )sin 0^ 


(280) 


(281) 


where the gimbal acceleration terms are again neglected. 
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Necessary gimbal torques are developed from 


G x -^xSx I 

G y = _I y g y sec ^g 


(282) 


and require the measurement of the gimbal angle Tg. 
c. Reaction Jets 

The reaction jet system comprises the propellant, oxidizer, 
engines, and tankage weight necessary for the spacecraft's control. To 
arrive at weight estimates for such a system it was assumed that the 
usable specific impulse considering engine efficiency, expulsion 
efficiency, and ullage would be 290 lbfsec/lbm and that the propellants 
would be storable hypergolics housed in tanks with positive feed expul- 
sion diaphragms. Manufacturer's data* can then be extrapolated (ref. 33 ) 
to yield the idealized total system weight 


w T = 0.0101 i t °* 9 12 (283) 

which is presented in figure 32. Here I T is the total impulse in lb sec. 
This impulse may be written in terms of the total momentum envelope for 
the spacecraft, giving 


i, . (Sac + M + 5b* 

\}x V 1 Z /T 


( 284 ) 


where Z x , Zy, and Z z are the moment arms about the X, Y, and Z body 
axes. 


*Reproduced and used by courtesy of the Minneapolis-Honeywell 
Regulator Company. 
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Figure 32.- Variation of total reaction jet system weight 
with required total impulse. 

An assessment of the weight of the reaction jet system thus 
requires the development of total momentum envelopes for the spacecraft 
mission. The momentum for particular disturbances must be determined 
from (25l) and the resulting momentum components along the spacecraft's 
axes must then be multiplied by the probable frequency of occurrence of 
each disturbance. By repeating this process for all disturbances and 
summing the individual momenta along each axis, a total momentum envelope 
per sampling period is obtained. The weight crossover time between 
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momentum storage and reaction control system can now be established from 
figure 32 , since the total momentum envelope per unit time has been 
developed. If the mission time exceeds this crossover time, momentum 
storage systems should be selected for damping of the spacecraft's 

motion. 

If reaction- jet systems are chosen, they can be combined with a 
mass-balancing system which compensates for any constant products of 
inertia resulting from crew motion or cargo transfer and eliminates limit 
cycling of the jets about the new principal moments of inertia. Such a 
system could, for example, pump the propellant to different positions 
within the spacecraft to obtain its control torques. Since the design 
of this mass-balancing system is very much dependent on the spacecraft 
geometry, it will not be considered here. 

In most cases, however, the control system will consist of both 
reaction jets and momentum storage systems. The jets then provide for 
attitude control and orbit keeping and the momentum storage system is 
used to damp any oscillatory motion of the spacecraft. Attitude control 
commands are now used to actuate the reaction jets while rate and rate 
integral commands provide control laws for the momentum storage system. 
The development of such combined systems is again dependent on specific 
spacecraft and disturbances and will not be attempted in this analysis. 

F. Comparison of Exact and Approximate Solution 
Two possible manned spacecraft were considered for a comparison of 
the results of the numerical integration of the exact equations of motion 
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and the results of the present analytical solution. These spacecraft 
were a cylindrical manned orbital research laboratory and a large 
hexagonal space station. 

1. Manned Orbital Research Laboratory - MORL 

The MORL is proposed as an earth-orbital laboratory in which scien- 
tific and engineering experiments could be conducted over extended time 
periods. The basic laboratory is designed to support a crew of six 
astronauts in a 200- nautical-mile orbit for up to 5 years. During 
spinning operation, the laboratory module and the last stage of its 
Saturn booster would remain attached by a system of cables and would 
rotate about a common mass center. The resultant centrifugal force 
would produce an effective gravity field in the manned module. A sketch 
of the corresponding MORL configuration is shown in figure 33- 



Figure 33.- Artist's sketch of possible manned orbital 
research laboratory. 
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Assumed inertia and mass characteristics (ref. 54) for this configuration 
are given in table 6. The inertia distribution is near-cylindrical, and 
the manned module and booster counterweight rotate about the Z axis at 
0.4 rad/sec. A gravitational acceleration, equal to one-fourth that at 
the earth's surface, acts on the manned module due to this rotational 
rate. 

The effects of various disturbances on MORL are summarized in table 7. 
The disturbances include residual rate and attitude errors after spinup, 
moments applied by an attitude control system valve failure in the open 
position, and several "worst— case” crew motions. These motions comprise 
step translation to an extreme position within the laboratory and linear 
oscillations which could result from trampoline exercise, ladder climbing, 
or floor pacing. A linear velocity of 4 ft/sec is selected for all 
oscillatory motions, and the entire crew of six is taken as a single 
equivalent mass with a mass factor Q of 3 6 slugs. 

Equations defining these disturbances are listed in the second 
column of the table. The resulting error l imi ts have been found by hand 
calculations of the analytical upper bounds and by extrapolation of the 
error time histories obtained from numerical integration of the exact 
equations of motion on an IBM 7094 computer. 

Both rate and attitude error limits are given in the table. The 
rate error limits range from about 0.004 rad/sec for the step inertia 
product to about 0.025 rad/sec for the step torques, and the attitude 
error limits vary from approximately 0.01 rad/sec for the step inertia 
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product to approximately 0.l8 rad/sec for the residual errors after 
spinup. Significantly, the errors caused by the periodic mass motions 
are several times greater than those produced by instantaneous motion 
to a final position. Predictions of maximum spacecraft errors due to 
crew motion must thus consider any periodic crew motions that may occur. 
Onboard experiments, which require high-accuracy control of the space- 
craft, may be adversely affected by the oscillatory crew motions and may 
require restriction of these motions. 

The approximate error limits developed from the analytical solutions 
show reasonable agreement with the maximum errors determined from the 
exact solutions. The deviations of the approximate error limits from 
the exact error limits are generally less than 20 percent of the exact 
error limits. These upper bounds of the spacecraft errors will thus 
give a conservative estimate of the effects of various disturbances 
and should be sufficient for initial engineering design applications. 

About 3 hours of hand calculations were required for the determina- 
tion of the approximate error limits, as compared with about 6 hours of 
data processing and computing for the calculation of the exact error 
limits. When one further considers the complexities associated with the 
programing and numerical integration of the exact equations of motion and 
the nonavailability of an IBM 709^ computer to many scientists, the 
advantages of the analytical results are apparent. 

Two of the solutions described in table 7 have been selected for a 
comparison of the actual error histories and the error histories given 
by the analytical solution. The disturbances are the step product of 
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inertia and the vertical mass oscillation. Both the uncontrolled and 
the controlled spacecraft motions were considered. Solutions for the 
uncontrolled case were obtained from numerical integration of the exact 
equations of motion and from evaluation of the error relations developed 
in this analysis. Solutions for the controlled case were obtained from 
numerical integration of both the linearized and the exact equations of 
motion. All calculations were carried out on an IBM 7094 computer. 

Single-axis control commands, which apply torques about the X or 
minimum- inertia axis of MORL, will be most efficient. One may accord- 
ingly set the nondimens ional rate damping gain k ly in ( 227 ) equal to 
zero. This yields 

k Ix = " T§- = - 2r D ( 285 ) 

At D 

as the nonaimens ional rate damping gain for the X axis. The corresponding 
rate integral gain, when used, is arbitrarily selected as 



to fall in the stable region of figure 27(d). 

For MORL, the time constant tp will be equated to one spin cycle 
or 20 sec. Values of the control constants are then 

Kq* = -15,300 

and 

K 2 x = ~ 2 >°^ 0 

The rate damping ratio becomes 

r D = O * 1 


L 
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and. state-of-the-art gyros (ref. 35) in the 500 to 1000 ft-lb-sec class 
(see fig. 3l) and jet hardware can provide the necessary control torques. 

The MORL response to the step product of inertia is illustrated in 
figures 3^ and 35* The uncontrolled rate and attitude errors are given 
in figure 3^« 



Figure 3^ •- Uncontrolled laboratory error histories for 
step inertia products. 


It is apparent that the exact and approximate solutions are virtually 
identical. Both the rate and the attitude errors are biased sinusoids. 
As expected from the trace analysis, the largest rate error occurs about 
the Y axis and the largest attitude error corresponds to rotation about 
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the X axis (see figs. 9 and 10). To the crew, the rate error appears 
as a minute rolling motion of the laboratory floor with a maximum 
amplitude of 0.6°. 

The controlled response of the laboratory to this disturbance is 
shown in figure 35* 




t. sec 
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Figure 35 -- Laboratory error histories for step inertia products 

and pure rate control. 

This figure, which corresponds to pure rate control about the minimum 
inertia axis, again gives the same results for both the exact and the 
approximate solution. The laboratory oscillation is reduced to steady 
coning in about three spin cycles. The constant residual rate errors 
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produce constant control torques which counteract the mass unbalance 
torques produced by the products of inertia. The effective disturbance 
torque is thus less than it was for the uncontrolled case, and the 
residual rate errors are correspondingly somewhat smaller than the con- 
stant components of the uncontrolled rate errors. The oscillatory terms 
in the uncontrolled rate errors are due to the acceleration terms asso- 
ciated with the introduction of the inertia products and tend to zero 
in the controlled error histories. 

The attitude errors for the damped rate errors become 


Op = q> r + i9 r = a o e“ i0t - — 232L ( 286 ) 

from (215). For the step products of inertia, the attitude errors result 
frcm changes in the body rates which do not affect the total spacecraft 
momentum. The contribution Oq of the transient oscillatory terms in 
the body rates to the attitude errors will tend to zero, and residual 
attitude errors are given by 


= 


^yr 

a 


and 


9 ^ = - 


foiCT 

a 


J 


(287) 


where and Hy r denote the residual body rates. Both the rate and 

attitude errors approach constant values for the theoretical solution. 

As predicted by the analysis of the controlled spacecraft character- 
istics, single-axis rate control is acceptable for noimal operation and 
experiments which do not require high-accuracy stabilization of the 
spacecraft . 


The MORL response to a vertical periodic motion of the entire crew 
is depicted in figures 36 and 37. Figure 36 illustrates the uncontrolled 
results. 


Approximate solution 
Exact solution 


Figure 36.- Uncontrolled laboratory error histories for 
vertical mass oscillation. 


The exact and approximate solutions check very closely. The rate and 
attitude errors now comprise a low-frequency, large-amplitude sinusoidal 
oscillation due to precession within the outer (A) ellipses and high- 
frequency, small-amplitude oscillations due to precession within the 
inner (p) ellipses (see figs. 12 and 13). The maximum errors are two 
to three times as high as the corresponding errors for the step inertia 





product. The laboratory floor also undergoes irregular rolling motions 
with maximum amplitudes of about 2°. Since the distance (ref. 34) from 
the center of rotation to the laboratory floor is approximately 50 feet 
this roll can produce a 2-foot total translation of the station floor 
and could present some difficulties to a moving astronaut within the 
laboratory. 

The controlled laboratory motion is presented in figure 37. 


- Approximate solution 

- Exact solution 



K Jx « -13,300 
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Figure 37.- Laboratory error histories for vertical mass oscillation 
and rate plus rate integral control. 


Control torques are applied about a single axis, but a combined rate and 
rate integral control law is used. The laboratory motion is quite similar 




l6o - 


to the uncontrolled motion, but exhibits damping of the free vibrations, 
as is apparent from the gradual decrease in the corresponding error terms. 
Since the periodic forcing function terms predominate, this type of 
response is to be expected. Agreement between the approximate and exact 
solutions is very good, and the small differences in the error histories 
can only be detected for the angle <p. 

The uncontrolled and controlled error histories developed from the 
approximate solution for the nonsymmetric MORL are practically coincident 
with the exact error histories for all the disturbances that have been 
examined. The analytical solution is thus a useful tool for the study 
of the nonsymmetric laboratory motion. 

2 . Large Manned Space Station 

A second possible type of manned rotating spacecraft is the large 
spinning space station, such as the 150- foot station which will be 
considered here. This station (ref. l), shown in figure 38 , has six 
cylindrical outer modules arranged in the shape of a hexagon. The outer 
modules are connected to a central hub and docking port by three spokes. 
Rotation about the maximum inertia axis provides artificial gravity for 
the living modules. The crew of this space station would vary from 
6 to 21 astronauts. 

Assumed characteristics for the 150-foot space station are listed 
in table 8. The inertia distribution approaches that of a flat disk 
and the spin rate is 3 rpm. The crew is taken as six astronauts with 
an effective mass factor Q of about 3 & slugs. 
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Figure 38 .- Artist’s sketch of possible 150- foot manned 

space station. 

Disturbance effects on this station are summarized in table 9* The 
disturbances are similar to the MORL disturbances, and periodic crew 
motions are simulated by motion of a single equivalent mass with a 
linear velocity of b ft/sec. The rate error limits range from approxi- 
mately 0.006 rad/sec for the step inertia products to approximately 
0.015 rad/sec for the circumferential mass motion. The attitude error 
limits vary from about 0.017 rad/sec for the step inertia products to 
about 0 .l 8 rad/sec for the residual errors. The errors due to periodic 
motions are considerably greater than those introduced by the step 


products of inertia. 
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Circumferential crew motions and the residual errors were chosen for 
a further comparison of the approximate and exact solutions. Both 
uncontrolled and controlled solutions were developed. 

Control torques were now applied about both station axes and the 
corresponding control gains were assumed to be equal. Referring to (226) 
and (227) , one notes that this yields 




as the nondimensional damping gains. The corresponding rate integral 
and attitude gains, when used, are selected as 

k 2x = “ k 2y = — 

and 

k 3x = " k 3y = k lx 

from the stable regions of figure 27(d) and 27(e). The physical control 
gains become 



Km = K2x = -222,817 

K 2x = " K 2y = -75,000 

and 

Kj x = -K 5y = -70,000 

from (210). A time constant of about three spin cycles or sec was 
selected to give the damping ratio 


r^ — 0.02 


Higher values of damping would require exorbitant control moment gyro 
and reaction jet control systems. Even the selected value will require 
gyros in the 5,000 to 10,000 ft-lb-sec class and will exceed the present 
state-of-the-art in gyro hardware (see fig. 3l) . Rapid jet damping, 
although feasible, will result in large fuel consumption. 

The advantages of single-axis control for nonsymmetric vehicles 
become obvious when one notes that the MORL, with about one-half the 
spin momentum of the 150-foot station, requires a control system that is 
an order of magnitude smaller. In addition, the MORL is able to achieve 
lower time constants and considerably better damping ratios. These 
results lead to the conclusion that nonsymmetric spacecraft, spinning 
about a maximum inertia axis, are preferable from the control standpoint 
and that single-axis stabilization about the minimum inertia axis can 
result in major control system weight savings for these spacecraft. 

The 150-foot station motion for the circumferential mass transfer 
is given in figure 39> The approximate and exact solutions are in good 
agreement, and the time histories exhibit slow oscillations. These 
oscillations (see figs. 21 and 22) consist of a large-amplitude sinusoid 
with the mass motion frequency p and small- amplitude sinusoids with, 
approximately, the precession frequency A. The angular deviation of 
the gravity vector has an amplitude of about 5° and appears as a cor- 
responding slow rolling motion to the crew. For the 150-foot space 
station this rolling motion produces a 4-foot oscillation of the station 


floor. 


Approximate solution 
Exact solution 


Figure 39.- Uncontrolled, station error histories for 
circumferential mass motion. 


The controlled station response is illustrated in figure 4o. Control 
is derived from rate plus rate integral commands, and the station rates 
are damped to a purely sinusoidal trace in approximately nine spin cycles 
The constant rate tern in the damped trace will eventually disappear 
under the action of the rate integral commands. Residual rate errors 
may then he expressed as 


ft xvr " _n r e 
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vhere denotes the half-amplitude of the residual rate. Frcm (160), 

this residual rate amplitude term is 

n Irzfa + p ) 2 

KX-p) 


so that 


Vr = " 


Irz( a + p) 


I(X - p) 


,ipt 


The corresponding residual rate error becomes 


Op = i 


^ p) 


I(\ - p) 


ipt 


and 


cp r = [ Jk_ ] 7 S - 

\<7 + p / 


= / \ ~ 
\a + p/ 



sin pt 


cos pt 


The residual terms correspond to the coefficients of e ipt in the uncon- 
trolled solution functions. The control system thus has to have little 
effect on the magnitude of the errors directly due to the constant 
circumferential mass motion as would be expected for the selected low 
value of the damping ratio. The approximate solution compares favorably 
with the exact solution for this example. 


1 66 - 






t. sec t, sec 

Figure 4o.- Station error histories for circumferential 
mass motion and rate plus rate integral control. 

Figure 4l presents the station response for residual rate and atti- 
tude errors. As anticipated from figure 15, the rate errors are simple 
sine and cosine curves. The attitude errors, following (133), are 
somewhat more complex sinusoids. The exact and approximate solutions 
checked to within three significant figures. 


Approximate and 
exact solution 
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Figure 4l.- Uncontrolled station error histories for 

residual errors. 


The controlled station motion with the residual errors is shown in 
figure 42. Rate plus attitude control commands are now employed, and 
the spacecraft completes the required 10° reorientation about two axes 
in approximately l6 spin cycles. The analytical and exact solutions 
again were identical. 
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Figure b-2.- Station error histories for residual errors 
and rate plus attitude control. 


In summarizing the comparison, one may conclude that the analytical 
solution was in excellent agreement with the exact solution for all cases 
considered. Since the spacecraft used in the comparison are typical 
examples of future rotating manned spacecraft, the analytical solution 
should he valid for the determination of the dynamics and control of most 
such spacecraft. Analytical results for unmanned spacecraft, which may 
have larger torque disturbances and residual errors but have few or no 
inertia changes, should also be acceptable. The time history data 
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obtained for the step torques and residual errors were accurate to three 
places and an increase in these disturbances should not appreciably 
degrade the results in the linear range. The analytical solution thus 
offers a simpler, more economical, and more direct means of assessing 
the effects of various disturbances and spacecraft characteristics on 
the spacecraft motion than the conputer runs. The insight into the 
mechanics of motion, that is gained from the error formulation developed 
in this analysis, should be of major value to future -work on the dyna mi cs 
of arbitrary rotating spacecraft. 
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IX. ANAIYSI8 OP SPINUP AND DESPIN MODE 


A. Governing Equations 

For the present application of the governing equations for the 
spinup and despin mode, the spacecraft disturbances are assumed to be 
restricted so that no internal mass movements occur and no moments are 
exerted about the spacecraft X and Y axes. This yields 


and 



(290) 


* 



Mz dt 


dt 


(291) 


from (8) and (l2). 

Spinup and despin moments about the Z axis will be assumed to be 
provided by constant- thrust, pulse-modulated jets (refs. 1, 34, 36). 
Since the control of the spacecraft during this mode is quite straight- 
forward, the main problem is the selection of a spinup and despin 
technique which minimizes the associated fuel consumption for rigid 
and extensible spacecraft. 


B. Rigid Spacecraft 

For the rigid spacecraft configurations, such as the large hexagonal 
space station, the jet moment arms remain constant. Spinup and despin 
fuel is thus given by 

W = - -^zf^zf 

SU z zf*SP z zfISP 


( 292 ) 


- 171 - 


where the subscript f denotes conditions after completion of the spinup 
maneuver. The simplest spinup technique would apply continuous thrust or 
constant-width thrust pulses until the desired spin rate is reached. 

C. Extensible Spacecraft 

1. Mathematical Model 

For extensible spacecraft configurations, such as the MORL, the fuel 
calculation becomes somewhat more difficult. As an example, consider 
the sketch below. 



Here the spacecraft consists of a manned module with mass and a 

counterweight module with mass m c . The two modules are connected by 
a flexible cable or strut arrangement, which is extended to produce a 
large rotational radius. The distance between the module mass centers 
is designated l j the offset distance between the thrust P m and the 
manned module mass center is and the offset distance between the 

thrust P c and the counterweight module mass center is Z c . 


These 


- 172 - 


offset distances yield a larger moment arm and may be required to main- 
tain the spinup thrust line normal to the line connecting the mass 
centers. To minimize oscillations of the manned module about its Z axis, 
some type of rate damping should also be provided. Rate damping moments 
can be supplied by a small reaction wheel or passive dampers. The 
individual module oscillations about their respective mass centers will, 
however, be neglected for the spinup fuel calculations. 

Spinup thrusts may be produced by jets on the manned module, by jets 
on the counterweight module, or by jets on both modules. The first method 
is preferable when the counterweight module mass exceeds the manned module 
mass, and the second method is preferable when the manned module mass 
exceeds the counterweight module mass. The third method may be used if 
jets are mounted on both the manned and counterweight module ; a pure 
couple about the spacecraft mass center can now be produced by selecting 



and 



( 293 ) 


Spinup thrusts would be simultaneously applied to both modules for this 
method. 

All three of the thrusting methods lead to effective moment and 
moment arm relations of the form 
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and 

“z = °3 Z z 



where 


and 


« *vPc + Pc% 

JJ j = 

( p m + p c)(®m + m c) 

D = p m^m + p c^c \ 

2 Pm + Pc 

% ’ p m + p o 


for the spacecraft model. 

The spin inertia I z can be written as 



(29^) 


(295) 


I z = I 0 - + Dcji 



an d I m , I c denote the respective module inertias referred to the 


module mass centers. 


2 .^ C. 
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The above equations now allow the simple formulation of the total 
fuel consumptions for different spinup techniques. 

2. Spinup and Despin Techniques 

Spinup will be assumed to occur in the following manner. While 
rigidly coupled, the two modules are brought to an ang ular rate Q zi . 

The modules are then separated by extending the flexible module connector 
under action of the centrifugal force. During this extension process 
the spacecraft momentum, spin-rate, or spinup thrust may be held constant. 
After the full cable extension is reached, the two modules are spun up 
to the final spacecraft spin speed. Despin will require this sequence 
in reverse order. 

The fuel required for spinup or despin can be expressed as 



where the subscript i denotes conditions after the initial spinup to 
^zi? ‘the subscript e denotes conditions after the extension; and the 
subscript f again denotes final conditions. 

Three characteristic spinup techniques will be considered here. 
These involve extension with constant momentum ^ z ±^ z ±> constant spin 
rate and continuous thrust. 
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a. Constant-Momentum Extension 

Perhaps the simplest spinup technique is one where no spinup 
thrust is applied during the extension. The modules are allowed to 
separate as desired while the spin speed automatically decreases to 
maintain the angular momentum constant. After achievement of the 
desired extension the spinup jets are again actuated. 

The fuel consumption now becomes 


W, 


SU 


= ikdke 


l zt J SP 


1 + 



( 299 ) 


from (298). The rate of cable (or strut) extension does not effect this 
fuel consumption and may be varied arbitrarily to maintain the cables 
in tension during the extension. The fuel consumption is minimized by 
selecting the smallest value of which will yield sufficient cable 

tension at the completion of the extension, 
b. Constant-Rate Extension 

For this spinup technique, the spacecraft spin speed is maintained 
at its initial value 2 z j_ throughout the extension. The modules are 
again allowed to separate until the final extension is reached, and the 
spacecraft is then brought to its final spin speed 

The fuel consumptions relation (298) reduces to 


W, 


SU 


I zf fl zf 
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^zi) 


( 300 ) 
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for this case. Fuel consumption now is minimized by selecting the 
smallest value of which will yield sufficient cable tension at 

the beginning of the extension. 

The spinup thrust for the constant-rate extension is established 
by the requirement that the rate of change of angular velocity due to 
operation of the thruster must be greater than that due to the rate of 
extension or retraction. When the extension rate is maintained at a 
constant value 1„ 9 this condition can be expressed as 


D5 ^ i ^ z I ^^5 ~ 


(301) 


The minimnm thrust is thus directly dependent on the product of the 
extension rate and the initial spin rate. Corresponding , conditions for 
variable extension rates may be developed from (290), if the time varia- 
tion of l z is known. 

c. Continuous-Thrust Extension 

Another possible spinup technique would involve continuous 
thrusting during the spinup. This brute-force technique will require 
rapid extension of the oable modules to be efficient, but will be simpler 
to implement than the constant-rate extension. 

The associated fuel consumption is given by 


*80 - ^ ^ ^ 
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( 302 ) 
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To evaluate (302), a time history for either 0. z or l z during the 
extension must he selected. Values of ft ze and t e can then be 
developed from this time history and ( 290 ), if one recalls that 
Z z = Z zf and - ^z = fi ze when t = t e . 

If the extension rate l z is constant, then 


"\ 


%e _ ~ — 
I zf 


• 

Izi^zi + I^^ziC'te - li) + — (te “ ti) 


and 
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Substitution of (303) into (302) now yields 
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and fuel consumption is optimized by selecting the smallest value of 
and the largest value of l z which will avoid cable slacking 
during the extension. 

d. Comparison of Extension Techniques 

To compare the different extension techniques, note that the fuel 
consumption for each technique is expressed as the ideal fuel consumption 
at full extension plus an incremental fuel consumption for the extension 
process. From the ratio of these incremental fuel consumptions for (299) 
and (300 ) , it follows that 
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is the condition corresponding to equality of the constant-rate and 
constant-momentum fuel consumption. Similarly the relation 
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must hold for equality of the constant-rate and const ant -momentum fuel 
consumption. 

The initial spin rates (fi zi ) CT and (^ z i) CR for the constant- 
thrust and constant-rate extensions should produce equal centrifugal 
forces to start the extension and are both assumed to be equal to the 

value (^zi)cR,CT' ■^ le spin rate (^zi)cM ^ or constant-momentum 

extension is greater than or equal to gip* 

The equations (305) and (506) are represented graphically in 
figure 44. This figure allows the direct selection of the most economical 
spinup technique for a particular spacecraft as a function of an extension 

1 rr -? 

length ratio , a moment arm ratio , a thrust ratio 

*zf *zi 
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, and a momentum ratio 
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To 


make use of the figure, one first selects the extension parameters 

o 

Dp, D^, Z 2 , (Dzi)cR,CT> and (^i)^* The selected parameters and 

the initial spacecraft characteristics determine values for the moment 
arm ratio, the thrust ratio, and the momentum ratio. 
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Figure 14.- Fuel criteria for extension technique selection. 

For a particular value of the extension length ratio, the moment 
arm ratio locates a fuel criteria point corresponding to the constant- 
rate extension and the thrust ratio locates a fuel criteria point cor- 
responding to the continuous-thrust extension. The momentum ratio gives 
a third fuel criteria point for the const ant -momentum extension on the 
ordinate. The lowest of the three fuel criteria values indicates the 
technique which will yield the lowest fuel consumption. 
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The special case, where the minimum centrifugal force for all 
three techniques is equal, is of interest since the minimum cable tension 
will usually determine the extension parameters. For this case, one has 




for the continuous-thrust extension. 

From (307) and (308), one may write 


Constant-momentum, fuel criteria + l zf 
Continuous- thrust fuel criteria 2 z i 

and the continuous-thrust extension is now always more economical than 

the constant-momentum extension. Only the constant-rate extension and 

the continuous- thrust extension need to be compared for this case. 

As an example, consider the Manned Orbital Research Laboratory 

described in table 6. The assumed extension parameters for this 



spacecraft are 


- l8l - 


l ± = 37.5 ft 
If = 137-5 ft 
(^zi)cR,CT = O* 1 rad/sec 
fl z f =0.4 rad/sec 
Dj = P m = 100 lb 
D 2 = l m = 4 ft 

and a spinup technique vhich will yield the lowest fuel consumption is 
desired. The cable tension must be greater than or equal to its initial 
value during the extension and the spinup time is immaterial. 

Since the minimum centrifugal force must be equal for all three 
techniques, the constant-momentum extension may be disregarded. From 
the given extension parameters and (30l) and (308), one has 

( ^z) CT = 0m2 5 2 ft/sec 

and 

( Z z )cR = ^*^-72 ft /sec 

as the respective maximum extension rates for the continuous-thrust and 
constant-rate spinups. 

The extension ratios now are 

M = 0.354 
Z zf 

= 0.254 
*zi 
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and 




^(^ct^^CRjCT 


1.02 


By referring to figure 44, one notes that 

(Fuel criteria ) =1.0 

(Fuel criteria)^ = 1.6 

Since the fuel criterion for the continuous-thrust extension is con- 
siderably lower than that for the constant-rate extension, it follows 
that the continuous-thrust spinup technique will require the least 
spinup fuel for this example. 

The actual fuel consumption values for the example, as computed 
from (299), (300), and (304) for a specific impulse of 290 lbf sec/lbm, 
are 


W CT = 245 1b 
w CR = 256 lb 
Wcm “ 276 lb 

The continuous-thrust spinup requires approximately 12 percent more fuel 
than the ideal spinup at full extension and the constant-momentum and 
constant-rate spinups require approximately 26 percent and 17 percent 
more fuel than the ideal value of 219 pounds. A saving of about 22 pounds 
of fuel can thus be realized for each spinup and despin cycle by selection 
of an optimum spinup technique for this example. 
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If desired, the response of the spacecraft to internal mass 
movements and external disturbance torques may be included in the 
analysis of the spinup and despin mode by using (8) and (l2) as the 
governing equations of motion. Supplementary linearized equations of 
motions for the relative module oscillations can be incorporated in such 


an analysis. 



X. CONCIZJSIOHS 


An approximate solution of the equations of motion of arbitrary- 
rotating spacecraft with variable disturbance functions has been developed 
on the basis of small changes in the spacecraft body rates, Euler angles, 
and inertia terms. Complex representations have been used to define 
spacecraft and rate errors induced by the disturbance functions, and the 
solutions for the time history components and total error vectors have 
been examined for both uncontrolled and controlled spacecraft. 

The results of this analysis have led to the following conclusions: 

A. A comparison of the present analytical solution and solutions 
obtained by numerical integration of the exact equations of motion for 
two typical manned spacecraft has shown that the analytical solution is 
in excellent agreement with the exact solution for the small angle and 
rate regime. The analytical solution provides a simpler, more economical, 
and more direct method of assessing the effects of various disturbances 
and spacecraft characteristics on the spacecraft motion and allows an 
insight into the mechanics of motion which cannot be derived from the 
numerical solution. 

B. Analytical upper limits of the rate and attitude errors induced 
by various disturbances are in reasonable agreement with the maximum 
errors found by interpolation of the numerical data. These upper limits 
should suffice for first estimates of the effect of the disturbances 

on the spacecraft motion. 

C. The spacecraft inertia distribution was found to have a 
significant effect on the spacecraft motion for equal disturbance 


- 185 - 


characteristics. Spacecraft, whose inertia distribution approached that 
of a flat disk, exhibit considerably more inherent stability than slender, 
near-cylindrical spacecraft spinning about a maximum or minimum axis of 
inertia. However, for most practical cases the disturbance characteristics 
are directly related to the inertia distribution, so that the error bounds 
for spacecraft with different inertia distributions will tend to be 
similar. 

D. Periodic mass motions within the spacecraft may result in rate 
and attitude errors, which are several times greater than those predicted 
for worst-case step products of inertia. For equal disturbance charac- 
teristics, the largest errors resulted frcm circumferential mass motion 
in the direction of spin. Respectively, smaller errors were produced 

by radial mass oscillations in an offset spin plane and vertical mass 
oscillations parallel to the spin axis. Motions of the crew such as 
trampoline exercise, ladder climbing, or periodic translations along 
the spacecraft floor should be carefully examined to dete rmin e their 
impact on the spacecraft motion. 

E. The spacecraft errors indicated instability trends, when the 
spin axis became an intermediate axis of inertia during a mass motion 

and when the periodic motions took place with the precession frequency A. 
Mass motions falling in these two categories should be avoided. 

F. An investigation of possible control techniques revealed that 
pure rate control and rate plus rate integral control would provide 
adequate damping of the spacecraft errors induced by internal disturb- 
ances. Initial attitude errors and attitude errors induced by external 
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disturbances cannot be eliminated by these control techniques and will 
require rate plus attitude control. 

G. Single-axis control was found to be acceptable for all control 
techniques and allovs major reductions in the control gains and control 
system weight for near-cylindrical configurations spinning about a 
maximum axis of inertia. 

H. Optimization of the spinup and extension technique for cable- 
or strut-connected spacecraft modules can lead to appreciable fuel 
savings for the extension and retraction process. Comparison of con- 
tinuous thrust, constant rate, and const ant -momentum extensions for an 
example spacecraft indicated that 22 lbs or 5 percent of the ideal spinup 
and despin fuel could be saved by use of a continuous-thrust extension 


technique . 
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XIV. APPENDIX A 

DEVELOPMENT OF THE LINEARIZED EQUATIONS OF MOTION 

The rotating spacecraft will be considered as the system of parti- 
cles shown in figure 1. A set of X Y Z axes fixed to the spacecraft is 
used to describe the rotational motion of the spacecraft with respect 
to a set of Xj Yj Zj axes which translate without rotation in inertial 
space and which remain parallel to a set of Xp Yp Zp axes fixed in 
inertial space. The general moment equation (ref. 6) about the origin 
of the X Y Z coordinate system is then 

(A- 1 ) 

It will be assumed that the system mass does not change during the time 
periods of interest so that 

“ = X ** x (A_2) 

The absolute vector acceleration Rj is given by 

T? 7? 7? — » — » — » —> — > 

Rj = Rq + rj + fl X rj + X fj + n x (n X r^) (A-3) 

and substitution of (A-3) in (A-2) yields 



(A-M 
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M = m s r s X Rq + 


[ ft x r j ) + y r j X m j ( h x r j ) 


+ ) x®j(« x V + 


X " X 


x m j ( ft x r j ] 


V 7? /-* ^ \ V -» " 

+ j n x m j( r j x r j) + j r j x m j r j 


(A-5) 


The acceleration of the origin Rq is found from the general force 
(ref. 6) equation 


fr? 7 ^ 7 * — > -* t* _> 

P = m s R 0 + r s + ft X r s + 2fl X r s + X (fi X r s ) (A-6) 


and the first term of (A-5) may now be written as 


— > — > 

-* 77 -» 


3 r s X H 0 = r s X P - r s X mg(ft x ir s ) - 7 S x m s (i2 X r s ) 


r q x m g ft x (ft 


X?.)] - 


r s X m s (fi s x r s ) - r s X m g r g 


= r'g x P - r s x m g (ft x r s ) - r g X nig (ft x ? g ) 


- r’a X m g (n X r s ) - ft X r, X m jt 


ft x r s ) - 
— 


ft X m s (r g X r g ) 


“ x m g r g 


(A-7) 
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The equation of motion becomes 


— » -» —> 

M = To X P + 


V 7» -» 7* -a> -h> 

/ r j x X rj) - r s X m s (fi x r s ) 


V X mj(fi X r^) - r s X m g (fi X ? S )J + £ 


V -» 

Tj X nij(fi X fj) 


? s X m s (fi X r s ) 


— > 

+ fi X 


V -» -» 

r j X X r A - r s x m s( n x r s)| 


+ < n x 


V -1 


m j r j x r j - m s r s X r - I + 




r j X n ,] r j - r s x “s' 


(A-8) 


where the vector from the center of mass to the origin is 


r s 


- V ^ 

A m R 


j 


(A-9) 


To reduce (A-8) to a more useful form, the particle system will be 
represented as a large mass associated with the spacecraft and fixed 
with respect to the X Y Z axes and n smaller masses which move rela- 
tive to the X Y Z axes. The rigid-body angular momentum vectors of the 

spacecraft, the n moving masses, and the spacecraft mass center will 
— ^ — > — > 

be designated as H 0 , H n , and H s , respectively. The rigid-body angular 


momentum of the system H referred to the system center of mass is then 


H " Hq + ^ “ H s 


(A-10) 
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where 

— ) > — > 

H g = r g x m s (n s X r s ) 

n 

H n = 7 rj x x rj) 

5 ) = 7 Xm j ( ^ x 7 } 

j=n+l 

and (A-8) can be rewritten as 




dH J -» 

M = r X P + -i-i + /fix 
s dt 


nn 

“ 3 > T* 

/ m j r j X r j - “ s r s X- r s 

J=i 


i 

j=i 


r j X m j r j - r s x m s r s 


(A-ll) 


(A-12) 


with 


? . V 51 ?. 

s U m B r 0 


J=1 


(A- 13) 


In component form, one obtains 

M x ~ ys ? z ~ z s p y + * ^xy^y ~ ^xz 4 \ + " ^xy^y " ^xz^z 

- ^(iyi^y - Iyz^z " 3-yx^x) + ^y(^z^z “ ^zx^x " ^-zy^y) 


n 


I m j(x 3 y 3 - y 3 i 3 ) 
L)=i 

“n 

pi 

"*n 

7 m^yjZj-zjyj) 

M 


m s ( x sys - y s x s) 


m s( x s^s “ z s x s) 


m s (y s z s - z sy B ) 


Qy 


>n„ 


(A-l4) 
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^y z s^x “ x s P 2 + ly^y - lyz^z “ -^yx^x + ^y^V " ^yz^z “ ^yx^x 


J=1 


-^zx^z “ I zy fl y) + 

%(lx% “ Ixy% 

(yjij - x jjj) 

- 

m s (y s x s - x s y s ) 

(yj'j - z «> 

- 

m s(y s z s - z s y s ) 

(z j*j - 

- 

m s ( z s x s - x s z s) 


fi„ 


(A- 15) 


and 


M z x s^y “ y s^x + ^-z^z ~ -*-zx^x ” -^-zy^V + ^z^z ” -^zx^ " ^zy^ 


" ^y(-*-x^x " ■’•xy^V “ ^xz^z) + ^x(ly^ " ^yz^z " lyx^x) 



m s ( z s x s “ 

x s z s) 

- 


r- 

- 

m s ( z s y s - 

y s z s ) 

p 

- 

m s ( x sy s - 

y s x s ) 


fix 

fiy 




(A-16) 



- 198 - 


whore 


■^x ^xo + / m j(yj + z j ) * m s(y s + z s ) 

a 


x y = J yo + 7 m j( x j 2 + z j 2 ) - m s( x s 2 + z s 2 


m R (x R + z K ) 


J z = x zo + ) m j( x j 2 + yj 2 ) - m s( x s 2 + y s 2 ) 


l XZ = V m j( x jZj) " m s( x s z s) 

i=i 


J yz = /, " m s(y s z s) 


m j( x jyj) - m s( x s y s ) 


I x = 2 / + z j z j) - m s (y s y s + z s z s ) 


ly = 2\) + ZjZj) - m g (x s x s + z g z s ) 


I z = 2 ) + y^yj) - m s ( x s x s + y s y s ) 
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X m «5 (x «A + " m s( x s z s + z s x s) 

j=l 



X m j + - m s (y s z s + z E y s ) 

j=i 



X + y jV - m s( x s y s + y s x s) 

j=i 


f (A- 18 ) 


J 


The coordinates of the jth moving mass are xj, yj, Zj, and the 
coordinates of the mass center are 

“V 

0 

0=1 



7s = L % } 


n 



(A-19) 


For the special case when only one mass m with coordinates x, 
y, z is moving with respect to the spacecraft, the equations of motion 
yield 


m f * • • • • • • 

Mx = m^ - ( yE>z “ + I X°X " ^xy^y " ^z^z + *x^ " I xy J V “ I xz fi z 

~ ^z ( ly^V “ ^yz^z " lyx%) + ^(iz^z “ Izx^x “ Izy^V) 

+ Q[(xy - yx)fly + (xz - zx)^ + (yz - zy) 


(A-20) 
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My m s ( zP x " *Pz) + I y^y - lyz^z “ lyx^x + iy*V " iyz^z - 


and 


- ^x( I z fi z - I zx n x - x zy^y) + ^(^x^x “ ^xy^y - Ixz°z) 
+ Q^yx - xy)% + (yz - zy)j22 + (zx - xz)j 


Mz m s ^y ■ yP x) + I z^z ■ X zx^x “ I zy^y + ^z fl z _ 


i V^ I x^ " I xy fl y “ ^Z^z) + ^x^y^y - I; 


sl - 


y z c *z 


P ZX^X ■ 

lyx^x) 


+ Q 


with 


(zi - xzjft*. + (zy - yz)0y + (xy - yi)J 


T x = ^O + Q(y 2 + z 2 ) 


I y = -Vo + Q ( x2 + z2 ) 


and 


I z = x zo + Q(x 2 + y 2 ) 
Ixz = Q(xz) 
lyz = Q(y z ) 

P xy = QC 3 ^) 

ix = 2Q(yy + z£) 

Iy = 2Q(xx + zz) 

I z = 2Q(xx + yy) 




\ 


%x^x 

(A-21) 

p zy*V 

(A-22) 
( A-23) 


(A- 2b) 
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I X Z = Q(xz + zx) . (A-2h) 

iyz = Q(yz + zy) 

ixy = Q(xy + yx) j 


where Q is given by 


Q = 


m(m s - m) 


m c 


(a-25) 


The spacecraft equations of motion (A-l4) to (A-l6) can be solved 
for the body rates ^y, and ^z« The motion of the rotating space- 
craft is then defined in terms of the modified Euler angles \|r, 9, 

and cp. These angles, as shown in figure 2, relate the moving body 
axes X Y Z to the intermediate reference axes Xj Yj Zj. From the figure 
one notes that the time derivatives of the Euler angles are 


<P = ftjj. + fly tan 9 sin cp + ^ tan 9 cos cp > 

9 = fiy cos cp - flz Sin cp \ (A-26) 

^ cos cp sec 9 + sin cp sec 9 


The Euler angles found from (A-26) and the body rates found from (A-l4) 
to (A— l6 ) completely define the rotational motion of the spinning 
spacecraft. 

For a large number of practical applications, one is concerned with 
motions involving small oscillations of the spacecraft spin axis from an 
equilibrium reference position. If it is assumed that the spacecraft 
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spins about its Z axis, that the Z axis is initially an axis of maximum 
or minimum inertia, and that the Zj axis is selected as the inertial 
reference, then 


sin cp = tan 9=9 cos cp = 1 

sin 6 = tan 0 = 9 cos 0 = 1 


(A-27) 


and 


«x « 

Oy « A Z 


(A-28) 


for the small oscillation regime. Consistent with these assumptions, 
one can consider the variable inertia terms to be sufficiently small in 
comparison with the spacecraft moments of inertia so that they may be 
neglected when multiplied by the oscillatory body rates or any angular 
accelerations. 

The reduction of the nonlinear governing equations to linear approxi- 
mations can best be accomplished by first converting (A-l4) through (A-l6) 
to nondimens ional form. As was done in reference k, one may introduce a 
nondime ns ional time t, an inertia term ^pq> and a nondimens ional 
mass Pj by 


T = ft zo t 


e pq - 


x pq. 

^zO 


<3 m„ 


(A-29) 


where p and q. range over x, y, and z. The remaining nondimensional 


terms are then 
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ca_ = —2— 
^ «zo 


oL = - 2 E_ 

^ V 2 


IzO^o 2 


^zo V/^zO^s 


u j = x j It 



^ ^zol/^zO 


u < = JjL /^L. 
J n zoV Iz0 


v j = yj 



m s 

W* = z i / 


and the nondimens ion al equations of motion become 


^ “ v s p z ~ w s p y + e xf^x ~ G xy^y " ^xz^z “ € x a bc " e xy^y " € xz co z 
" a) z( € y a) y “ € yz CD z " € yxJ-°x) + ^^yC^zPz ~ ^zxf^bc " ^zy^^y) 


l 

j =i 


Y Hj(ujVj - VjUj) 

i=i 


- ( u s^s - v s^s)> “V 


- (U S W S - VgUg) 



^j (v Y’j ” Vj) " (v s^s - w s V s) 
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^y - w sPx “ u sPz + e j4 J y ~ c yz^z " ' £ yy^x. + £/°y - e-yzPz “ ^yx^bc 
- (^(e^Pz ~ e zx^x “ e zy x> y) + tu z( € xPbc " € xy aj y " ^xz^z) 



J=1 


- u jV 


- < v s“s - u s' , s>) <“x 


n 


V HjCvjVj - WjVj) 

J=i 


n 

) ^(yjUj - UjWj) 


<3=1 


- (v s w s - W s v s ) > 03 Z 


- ( w s^s - U s^s) 


(A-32) 


and 


^z “ u s p y ~ v s p x + e z^z " € zxf^x ~ e zy^y + G z a) z “ ~ e zy^y 


- “y( e A " e xy a V - G xz^z) + “bcC^V “ e yz^z " ^yx^) 



(v s u s - u s w s )> 



(A-33) 
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where 


€ x = c xo + 


n 




7 Hj(vj 2 ♦ „/) . (v/ + V S 2 ) 


r n 


e y e yo + 


7 “j (u / + W J 2) • (u s 2 + “s 2 > 

3=1 


€ z e zo 


n 


V Hj(uj 2 + Vj 2 ) - (u s 2 + V S 2 ) 
0=1 


\ (a-34) 


-xz 


n 


V Hj(ujWj) - (u s w s ) 


0=1 


-yz 


n 

v 

o=i 




- ( v s w s) 


'xy 


y njCujVj) - (u s v s ) 

<Fi 




In accordance with the small oscillation assumptions, take 9qp, 

* n2 2 

9 > 9 > Lp, Pp, and the variable inertia terms 

n n 

7 v i n i 2 ’ 7 

0=1 0=1 


^ v 3 
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v 

L 

0=1 


w 2 ' 


n 

7 

j-i 


“3 U 3 V 3' 


n 

V 1 

m 3 u 3 w 3» 

3=1 
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^0 V 0 W 0' 


0=1 


- 20 6 - 


and their derivatives to be of order A. Here A is restricted to be 
sufficiently small so that terms of higher order than A may be 
neglected in the governing equations. 

Thus, it follows that 



0[A) 



0 A 
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and 


^z ~i u s p y ~ v s p x, + . e z^°z ^zxf^x. “ ^zy 0 ^ + e z a) z " e zxf^x " e zy a) y 




0@ 00] 00 003 00] 00 003 003 


ax 


■y(^“ ^x/V, “ £cz“^) + 5y^z, - ^y^) 

003 003 o0] 003 00] 003 


n 


A * Vo* 




- ( w s^s - v s^s)/ “V 






+ 




< I ^Vj - Vi 1 

l Lj e ; 



(A-3T) 


00 
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or neglecting terms of higher order 


~ + ( € z " “ ^z^xz " ^yz 0 ^) 






I y ~ € y x> y ~ ( e z " e x) a bc? D z " ^zC^yz + ^z^z) 


+ 


r 


r n 



(v s w s 



J 


+ 






J 


L z = e A + e z co z 



(A-38) 


(A-39) 


(A-it0) 
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For the same range of disturbances, the Euler angles are given by 



and again neglecting terms of order higher than A 



(A- 4 l) 


(A- 42 ) 


If the spacecraft dynamics are well conditioned, the above equations 
should give reasonable results in the small angle and rate regime, for 
which 


-15° < 9 < 15° 
-15° < e < 15° 

-O.0685 < A < O.0685 
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Solutions for higher values of A will lead to correspondingly less 
accurate solutions. 

Since the form of the physical and nondimensional differential 
equations is identical, only the more convenient physical equations of 
motion will he used. These equations can then be written as 



Y“\Mx + %(ixz “ 


"" n 

y - x j 

x l 


1 

J= i 

- m s ( z s x s - X s z s) 

+ 

~ n 

v 
/ , 

- yj’ z V 



J=1 





- m s (z s y s - y s z s ) 


(A-43) 





Ixz%) 


+ 


m 



m s( z sys - ys z s) 

ftz + 

~n ^ 

y m j( x j z j ■ z j x j^ 

_ 


jFi 


and 





j=i 


(A-44) 


m s (y s x s 



(A-45) 
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where the moments of inertia I x and Iy, in general, may be approxi- 
mated by their initial values in (A-43) and (A-44), so that 


Ix - IxO 


x y - IyO + 


” -*-z0 


n 

7 mj(y Jo 2 + z io 2 ) 


j=l 


~ * s (y so 2 + z so 2 ) 


n 

7 m «j( x jo + z jo ) 


3=1 


' n 


m 


s( x so + z so ) / (A-46) 


7 m^Xj 2 + yj 2 ) 

3=1 


- m s (x s 2 + y s 2 ) 


J 


since the retention of variable inertia coefficients does not appear to 
add appreciably to the accuracy of the solutions. The other pertinent 
inertia terms are 


L xz 




~ n ~~ 



- m s( x s z s) 

j=i 


— n 


I m a^3 z j> 

- m s (y s z s) 

J=1 

) 


L yz 


and the required time derivatives of the inertia terms become 


L xz 


L yz 


Iz - 2 


(A-47) 
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m j( x j z j 

+ z j x j) 

- m s (x s z s + z s x s 


Lra. 


— 



~ n 





y 

m j(yj z 3 

+ z 3yj) 

- %(y s z s + z sys) 



- 



~ n 
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+ yjyj) 

- %( x s x s + y s y s ) 


j=i 




) 


(A-48) 


- 213 - 


The Euler angle relations are given by 

4 > - ^0 = 0 * 

0 + ft z cp = Qy. 

t = fiz 

If results should be required in nondimensional form, one need only 
apply the transformations (A-29) and (A-30) to the solutions. 
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XV. APPENDIX B 

PARTICULAR SOHJTION FUNCTIONS FOR THE 
UNCONTROLLED SPACECRAFT 


The forcing function for the differential equation with the 
stant moments of inertia approximation can he written as 


JL 

F = V AjFj(t) 

J=i 


and thus the particular solutions are 




r * A A J F o (t) 

j=l 


F = V AjFj(t) 


J=1 


F = 7 AjFj(t) 
j=l 


where the functions Fj(t), Fj(t), and F^(t) are given hy the 
Laplace transforms 


FjM 


mi 

'.<*>- mr m-g 

^■my-im- 34 


s) 


con- 


(B-l) 


(B-2) 


inverse 


(B-3) 
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A number of particular solution functions have been evaluated and are 
presented in tables 1-3. Other functions may be determined from 
(B-3) if needed. 

The unit step function and the unit impulse function, that occur 
in these tables, are defined as 



where t = 0 is the initial time of application of the step or impulse 
disturbance. 
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J 
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| i 
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i 
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TRI GONQMETRIC FUNCTIONS 



cos at - cos At 
































term l rz denotes the constant Inertia product 






























denotes the constant inertia product Qr 0 z ( 
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TABLE 6.- ASSUMED CHARACTERISTICS FOR MANNED ORBITAL 
RESEARCH LABORATORY 


Parameter 

Module values 

Total 

values 

Manned 

Counterweight 

2 

I x0 > slug-ft 

103,000 

30,000 

133,000 

IyO> Slug-ft 2 

90,500 

73,000 

7,393,412 

I Z Q, Slug-ft 2 

173,000 

73,000 

7,475,912 

m s , slugs 

1,220 

557 

1,777 

Q, slugs 

36 

— 

36 

cr, rad/sec 



0.4 
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TABLE 8.- ASSUMED CHARACTERISTICS FOR 150-FOOT SPACE STATION 


Parameter 

Total value 

2 

I x0 , slug-ft 

10,500,000 

I y0 , slug-ft 2 

10,500,000 

2 

I z0 , slug-ft 

15,000,000 

m s , slugs 

2,270 

Q, slugs 

36 

a, rad/sec 

0.31^ 






TABLE 9.- DISTURBANCE EFFECTS ON 150-FOOT SPACE STATION 
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